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PREFACE 


The  work  reported  herein  was  conducted  for  the  U.  S.  Army  Engineer 
District,  Walla  Walla  (NPW) , by  the  Structures  Laboratory  (SL)  of  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES).  Authorization 
for  this  investigation  was  given  in  DA  Form  2544,  dated  15  March  1978. 

The  material  property  information,  construction  history,  instru- 
mentation data,  and  the  Walla  Walla  version  of  the  temperature-calculation 
program  used  in  thi3  investigation  were  provided  to  WES  by  NPW  and  North 
Pacific  Division  Laboratory. 

This  investigation  was  performed  under  the  direction  of  Messrs.  B. 
Mather,  J.  M.  Scanlon,  and  J.  E.  McDonald,  SL.  Dr.  Tony  C.  Liu,  Messrs. 

R.  L.  Campbell,  and  A.  A.  Bombich  prepared  this  report. 

The  Commander  and  Director  of  WES  during  the  conduct  of  this  in- 
vestigation and  the  preparation  and  publication  of  this  report  was 
COL  J.  L.  Cannon,  CE.  Mr.  F.  R.  Brown  was  Technical  Director. 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiple 

By 

To  Obtain 

Btu  • inch  per  hour  • 
square  inch  * degree  F 

20.7688176 

watts  per  metre  • Kelvin 

Btu  per  hour  • square 
foot  • degree  Fahrenheit 

5.678263 

watts  per  square  metre  • 
Kelvin 

Btu  per  pound  (mass)  • 
degree  Fahrenheit 

4186.8 

joules  per  kilogram  • Kelvin 

calories  per  gram 

4.184 

kilojoules  per  kilogram 

cubic  yards 

0.7645549 

cubic  metres 

Fahrenheit  degrees 

5/9 

Celsius  degrees* 

feet 

0.3048 

metres 

inches 

0.0254 

metres 

inches  per  degree  Fahrenheit 

0.014111111 

metres  per  Kelvin 

pounds  (force)  per  square 
inch  per  minute 

114.91267 

pascals  per  second 

pounds  (force)  per  square 
inch  (psi) 

6894.757 

pascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

pounds  per  cubic  yard 

0.5932764 

kilograms  per  cubic  metre 

square  feet  per  hour 

0.0000258064 

square  metres  per  second 

' 


To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read- 
ings, use  the  following  formula:  C = (5/9) (F  - 32).  To  obtain  Kelvin 
(K)  readings,  use:  K * C + 273.15. 


i 
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VERIFICATION  OF  TEMPERATURE  AND  THERMAL 
STRESS  ANALYSIS  COMPUTER  PROGRAMS  FOR  MASS  CONCRETE  STRUCTURES 


PART  I : INTRODUCTION 


Background 

1.  Currently,  two  finite  element  computer  programs  are  being  used 
by  the  Corps  to  predict  temperature  distribution  and  resulting  thermal 
stresses  and  strains  in  mass  concrete  structures.  The  first  program,^"  de- 
veloped by  Dr.  E.  L.  Wilson  of  the  University  of  California  (UC)  at  Berkeley, 

calculates  temperatures  within  a mass  concrete  structure.  The  second 
2 

program,  written  by  Dr.  R.  S.  Sandhu.et  al. , also  at  UC  Berkeley,  calcu- 
lates the  thermal  stresses  and  strains  within  the  structure  involving 
incremental  construction  and  creep.  These  two  computer  programs  are  con- 
sidered to  be  the  most  effective  numerical  methods  for  calculating  tem- 
perature and  thermal  stresses/strains  in  the  mass  concrete  structures. 

Because  they  are  not  only  completely  general  with  respect  to  geometry,  mate- 
rial properties,  and  boundary  conditions,  but  a1 so  provide  the  capability 
of  simulating  incremental  construction  and  creep  of  mass  concrete  struc- 
tures, these  programs  have  been  used  for  concrete  temperature  control 
studies  for  several  Corps  mass  concrete  structures,  Including  Dworshak 
Dam.  However,  these  programs  have  never  been  verified  with  actual  field 
measured  data. 

2.  As  in  many  other  dams,  stress  and  strain  meters,  thermocouples, 
and  other  instruments  were  embedded  in  Dworshak  Dam  during  construction 
to  give  information  on  structural  behavior.  Measurements  have  been 
obtained  from  these  instruments  since  the  beginning  of  construction  in 
1968.  These  measured  data,  together  with  the  available  laboratory  in- 
formation on  concrete  thermal  and  mechanical  properties,  provide  a unique 
opportunity  for  verifying  the  computer  programs. 


4 


Purpose  and  Scope 


3.  The  purpose  of  this  investigation  was  to  verify  the  temperature 
and  thermal  stress/strain  computer  programs  by  comparing  the  analytical 
results  with  the  measured  data  obtained  from  the  embedded  gages  in  the 
Dworshak  Dam. 

4.  Because  of  the  funding  and  time  limitations,  the  analysis  of 
the  complete  dam,  taking  full  account  of  the  time-varying  thermal,  elas- 
tic, and  creep  properties  of  concrete,  was  not  considered  possible.  It 
was  decided  that  the  objective  of  this  investigation  could  be  achieved 
by  analysis  of  the  lower  region  of  the  Monolith  23,  including  foundation 
and  14  concrete  lifts,  during  the  period  between  26  August  1968  and 

15  February  1969. 

5.  History  of  construction,  weather  data,  and  available  material 
properties  of  Dworshak  Dam  were  used  for  computer  program  input.  Calcu- 
lated temperatures  and  stress/strain  history  were  compared  with  the 
measured  data.  The  validity  and  reliability  of  the  computer  programs 
are  discussed,  and  recommendations  for  further  development  of  tempera- 
ture and  thermal  stress  calculation  programs  are  presented. 


PART  II:  FINITE  ELEMENT  COMPUTER  PROGRAMS 


Temperature-Calculation  Program 

6.  This  program,  developed  by  Dr.  Wilson,  calculates  the  tempera- 
ture distribution  as  a function  of  time  within  a concrete  structure  as  It 
Is  being  constructed.  Each  lift  of  the  structure  may  be  placed  at  an 
arbitrary  time  and  temperature.  Insulating  forms  may  be  placed  or  removed 
from  the  surfaces.  The  external  air  temperature  and  temperature  of  the 
cooling  water  may  also  vary  with  time.  The  finite  element  technique 
coupled  with  a step-by-step  time  integration  procedure  is  used  as  the 
method  of  analysis.  Detailed  discussion  of  this  computer  program  is  given 
in  Reference  1. 

7.  Some  input  and  output  formats  of  Dr.  Wilson's  program  were 
modified  at  WES  for  the  convenience  of  the  users.  Both  the  WES  and  Walla 
Walla  versions  of  the  temperature-calculation  program  were  used  in  this 
investigation  and,  as  expected,  the  results  were  found  to  be  identical. 

The  WES  version  of  temperature  program  is  given  in  Appendix  A. 

Thermal  Stress/Strain  Calculation  Program 

8.  This  program,  developed  by  Dr.  Sandhu,  calculates  the  displace- 
ments at  each  node  and  the  strains  and  stresses  developed  in  each  element 
in  the  finite  element  model  due  to  thermal,  gravity,  and  other  external 
loads.  In  allowing  for  creep,  this  program  assumes  that  relaxation  of 
stress  takes  place  without  nodal  displacements  over  a small  time  incre- 
ment during  which  the  material  properties  do  not  change.  This  change  in 
stress  is  then  neutralized  by  releasing  the  constraints  and  treating  the 
stress  changes  as  residual  stresses.  Detailed  discussion  of  this  program 
is  given  in  Reference  2. 

2 

9.  In  Dr.  Sandhu' s original  program,  only  one  arbitrary  reference 
temperature  (stress-free  temperature)  was  specified  for  all  elements. 

The  stress-free  temperature  input  was  modified  in  the  WES  version 

of  the  thermal  stress  program  in  which  the  stress-free  temperature 
for  an  element  is  defined  as  the  temperature  at  8 hr  after  placement. 

6 


The  stress-free  temperatures  for  all  elements  are  calculated  in  the  WES 
version  of  the  temperature  calculation  program  and  stored  on  disc  file 
for  subsequent  input  to  the  thermal  stress  analysis. 

10.  The  WES  version  of  the  thermal  stress  analysis  program  was 
used  in  this  investigation  (Appendix  B). 
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PART  III:  DETERMINATION  OF  CONCRETE  TEMPERATURES 


Finite  Element  Model 

11.  The  finite  element  model  selected  for  this  inves -igation  is 

representative  of  a small  interior  portion  of  Monolith  23.  Figure  1 

shows  the  typical  cross  section  of  Monolith  23  and  Figure  2,  the  typical 

cooling  pipe  arrangement.  Due  to  considerations  of  symmetry,  only  a 
* 

30-in.  wide  slice  was  analyzed  as  shown  in  Figure  2. 

12.  The  height  of  the  finite  element  model  was  80  ft  (from  El  930 
to  El  1010),  representing  10  ft  of  foundation  and  14  lifts  of  concrete  at 
5 ft  per  each  lift.  The  plan  of  the  model  was  located  at  330  ft  from  the 
construction  base  line  and  center  of  the  monolith.  The  finite  element 
model,  consisting  of  576  nodes  and  475  elements,  is  shown  in  Figure  3. 

13.  Except  for  the  top  concrete  surfaces,  all  boundary  conditions 
were  assumed  to  be  adiabatic  because  negligible  boundary  heat  flow  was 
expected  in  the  interior  portion  of  a mass  concrete  structure. 

Material  Types 


14.  Four  adjustments  in  fly  ash  and  cement  contents  in  the  con- 
crete mixtures  were  made  between  August  and  December  of  1968.  Variations 
in  concrete  mixtures  would  affect  the  thermal  and  mechanical  properties 
of  concrete.  Therefore,  five  material  types  (one  foundation  material 
and  four  concrete  materials)  were  used  in  the  analysis.  A summary  of 
material  types  used  is  given  in  Table  1. 


A table  of  factors  for  converting  U.S.  customary  units  of  measurement 
to  metric  (SI)  units  is  given  on  page  3. 
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PIPES  PLACED  ON  TOP  OF  LIFT  PRIOR 
TO  PLACEMENT  OF  NEXT  LIFT 


OO LING  PIPE 


COOLING  PIPE 


TYPICAL  ELEMENT  OF  OAM 
FOR  TEMPERATURE  STRESS 


iOOLING  PIPE 


COOLING  PIPE 


COOLING  PIPE 


COOLING  PIPE 


•COOLING  PIPE 


COOLING  PIPE 


FOUNDATION 


Figure  2.  Typical  Arrangement  of  Cooling  Pipes 


Figure  3.  Finite  Element  Model 


Figure  3 (Continued) 
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Figure  3 (Continued) 
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3 (Continued) 


-a.3aaec  at 


Figure  3 (Continued) 
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Figure  3 (Continued) 


Figure  3 (Continued) 


Figure  3 (Concluded) 


Thermal  Properties  of  Concrete  and  Foundation 

15.  The  thermal  properties  of  concrete  and  foundation  required 

by  the  temperature  calculation  program  are  (a)  specific  heat,  (b)  density, 
(c)  thermal  conductivity,  and  (d)  adiabatic  temperature  rise  versus  age. 
The  values  of  thermal  properties  used  in  the  analysis  are  given  in  Tables 
2 and  3. 

16.  The  specific  heat,  density,  and  thermal  conductivity  of  mate- 
rial Type  2 were  obtained  from  the  test  results  of  mixture  No.  A-2  re- 
ported in  Appendix  C of  Reference  4.  The  thermal  properties  of  other 
concrete  material  types  used  in  the  analysis  were  interpolated  from  the 
available  test  data  given  in  Reference  4.  The  thermal  properties  of 
foundation  material  (material  Type  1)  were  assumed  to  be  identical  to 
material  Type  2. 

17.  The  adiabatic  temperature  rise  data  presented  in  Table  3 were 
derived  from  the  test  results  given  in  page  2.3-9  of  Reference  4.  The 
adiabatic  temperature  rise  data  beyond  28  days  were  obtained  by  extrapo- 
lation of  the  available  data  on  a semi-logarithmic  plot. 

Construction  Parameters 

Ambient  temperatures 

18.  The  mean  daily  ambient  temperature  at  Dworshak  Dam  from 
26  August  1968  to  15  February  1969,  during  which  the  construction  is 
being  simulated,  is  plotted  in  Figure  4.  This  information  was  provided 
by  the  Walla  Walla  District.  The  simplified  ambient  exposure  tempera- 
tures used  in  the  analysis  are  shown  in  Figure  5. 

Foundation  temperature 

19.  In  order  to  simulate  the  foundation  temperatures  in  the  model, 
the  temperature  of  the  entire  foundation  was  first  set  equal  to  65°F. 

It  was  assumed  that  the  temperature  at  the  bottom  of  the  foundation  would 
remain  constant  throughout  the  analysis  period.  Prior  to  the  beginning 
of  construction  simulation,  the  surface  of  the  foundation  was  exposed  to 
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Figure  4.  Mean  Daily  Ambient  Temperatures,  Dworshak  Dam 


Figure  5 Ambient  Temperatures  Used  in  the  Analysis 


ambient  temperature  and  calculation  was  made  to  equilibrate  foundation 
temperatures  between  ambient  at  the  surface  and  the  constant  temperature 
at  the  bottom. 

Concrete  placement  schedule 

20.  The  concrete  placement  schedule  used  in  the  analysis  is  shown 
in  Table  1,  which  was  obtained  from  Ref  ence  3.  As  can  be  seen,  the 
placement  rate  varied  between  4 and  13  days  per  lift  with  all  lifts  being 
5-ft  height.  The  first  lift  (El  940)  was  placed  on  27  August  1968,  and 
the  last  lift  simulated  in  the  analysis  (El  1010)  on  17  December  1968. 
Concrete  placement  temperature 

21.  At  Dworshak  Dam,  concrete  placing  temperature  was  specified 
between  40  and  45°F,  measured  20  minutes  after  mixing."^  Average  con- 
crete placement  temperature  of  42.5°F  was  assumed  in  the  analysis. 
Concrete  post-cooling 

22.  Post-cooling  of  all  concrete  was  required  in  Monolith  23  from 
the  base  of  the  dam  to  El  1175. 

23.  In  the  finite  element  model,  the  cooling  pipes  were  located 
at  nodes  67,  103,  139,  175,  211,  247,  283,  319,  355,  391,  427,  463,  499, 
535,  and  571. 

24.  The  inlet  and  average  outlet  cooling  water  temperatures  were 
41  and  53°F,  respectively.  Since  the  direction  of  flow  was  automatically 
reversed  every  12  hr,  the  average  value  of  inlet  and  outlet  temperatures 
was  used  in  the  analysis  (i.e.,  47°F). 

25.  The  water  was  flowing  in  the  pipe  at  the  start  of  placement 
in  a lift  and  generally  flowed  continually  for  21  days.  The  cooling 
water  removal  schedule  used  in  the  analysis  is  given  in  Table  4. 

26.  The  heat  transfer  coefficient  was  calculated  for  pipe  cooling 
based  on  the  following  equation:^ 


log  - - 2 
°e  r 


where  H = Heat  transfer  coefficient 

K * Conductivity  of  concrete 
= 0.113  Btu/hr-in.-°F 
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R - Mesh  size  In  the  vicinity  of  the  pipe 

■ 6 in. 

r - Radius  of  cooling  pipe 

■ 0.5  in. 

Because  of  symmetry  at  the  cooling  pipes,  the  mathematical  model  (Fig- 
ure 2)  included  only  concrete  elements  from  the  pipe  location  to  a dis- 
tance half-way  'to  the  next  cooling  pipe.  Therefore,  half  of  the 
calculated  H value  was  used  in  the  program  calculations. 

Surface  insulation 

27.  During  construction,  the  top  lift  surface  was  generally  covered 
with  insulating  mats  immediately  after  placing.  The  insulation  remained 
in  place  until  the  lift  surface  was  to  be  covered  with  concrete,  except 
for  the  time  during  which  the  joint  was  cleaned  with  high-pressure 

water  jet. 

28.  The  thermal  conductivities  of  insulation  materials  used  in  the 

3 

analysis  are  shown  in  Table  5. 

Wind  velocities 

29.  In  order  to  compute  the  convection  transfer  coefficient,  wind 
velocity  data  are  required.  The  wind  data  used  in  the  analysis  were  ob- 
tained from  Dworshak  Station  and  are  summarized  in  Table  6. 

Heat  transfer  coefficients 

30.  The  surface  heat  transfer  coefficients  used  in  the  analysis 
were  calculated  according  to  the  following  equation;^- 


where  H - Heat  tranfer  coefficient 

c = Conductance  of  insulation 

h'=  Convection  transfer  coefficient 
- 0.99  + 0.21V  for  V < 16  fps* 

V = wind  velocity 


From  ASHRAE  Handbook,  1977  Fundamentals. 
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31.  Based  on  the  data  given  in  Tables  5 and  6,  the  surface  heat 
transfer  coefficients  used  in  the  analysis  are  summarized  in  Table  7. 


Results  and  Discussion 


32.  The  analysis  was  performed  for  174  days,  from  26  August  1968 
to  15  February  1969.  Calculated  temperatures  were  output  at  one-day  in- 
tervals for  all  nodes  in  the  model  at  the  particular  stage  of  construc- 
tion. 

33.  In  general,  peak  temperatures  of  approximately  66°F  were  at- 
tained after  about  five  days.  Placing  of  precooled  overlying  lifts  re- 
sulted in  lower  temperatures  subsequent  to  their  placement.  After  at- 
tainment of  peak  temperature,  the  temperature  of  a lift  dropped  as  long 
as  cooling  was  continued.  When  cooling  of  a lift  was  stopped  at  21  days 
but  continued  on  the  overlying  lifts,  the  temperature  rose  gradually  to 
approximately  the  primary  peak  at  about  40  days,  and  this  temperature  re- 
mained essentially  constant  throughout  the  remaining  analytical  period. 
The  comparisons  between  calculated  and  measured  temperatures  are  made  at 
nodes  198,  342,  and  486,  where  thermocouples  2309B,  2319B,  and  2329B, 
respectively,  are  located.  As  can  be  seen  from  Figures  6,  7,  and  8, 
excellent  agreements  between  calculated  and  measured  temperatures  are 
obtained.  The  average  differences  between  calculated  and  measured  tem- 
peratures were  0.80°F,  1.56°F,  and  1.57°F  for  nodes  198,  342,  and  486, 
respectively.  These  small  discrepancies  were  due  probably  to  the  follow- 
ing reasons. 


a.  The  adiabatic  temperature  rise  data  used  in  the  analysis 
were  based  on  the  28-day  test  results . It  is  obvious  that 
heat  is  still  being  generated  for  an  extended  period  of 
time  after  the  28-day  adiabatic  temperature  rise  test 
period.  The  adiabatic  temperature  rise  data  beyond  28  days 
were  assumed  in  the  analysis  because  this  information  was 
not  available. 

b.  Average  concrete  placement  temperature  of  42.5°F  was  as- 
sumed in  the  analysis.  No  record  of  actual  concrete  tem- 
peratures at  placement  was  available. 

£.  The  average  cooling  water  temperature  was  used  in  the 

analysis  because  the  ectual  inlet  and  outlet  cooling  water 
temperatures  were  not  available. 


TINE  INTO  CONSTRUCTION  - OAVS 


<i.  The  thermocouple  2319B  clearly  indicated  that  the  concrete 
cooling  rate  was  too  rapid  at  node  342  , and  the  water 
flow  was  interrupted  between  Day  70  and  Day  75  (Figure  7). 
In  the  analysis,  however,  only  the  specified  cooling  flow 
was  used.  Discrepancies  between  measured  and  calculated 
temperatures  between  Day  70  and  75  were  to  be  expected. 

e.  Simplified  ambient  temperatures  and  wind  velocities  were 
used  in  the  analysis. 

f.  No  solar  radiation  was  considered. 


PART  IV:  DETERMINATION  OF  THERMAL  STRESSES 


Finite  Element  Model 

34.  The  same  finite  element  model  used  in  the  temperature  calcu- 
lation program  was  used.  Because  of  symmetry,  no  horizontal  displace- 
ment at  the  vertical  boundary  was  permitted  and  shear  force  was  assumed 
to  be  zero.  The  nodal  points  at  the  bottom  of  the  foundation  were  as- 
sumed fixed  in  all  directions. 

Properties  of  'Concrete  and  Foundation 

35.  The  properties  of  concrete  and  foundation  required  for  the 
thermal  stress  analysis  are  (a)  modulus  of  elasticity  versus  age,  (b) 
Poisson’s  ratio  versus  age,  (c)  coefficient  of  thermal  expansion,  (d) 
shear  foundation  factor,  and  (e)  creep  versus  age. 

36.  The  values  of  modulus  of  elasticity,  Poisson's  ratio,  and  co- 
efficient of  thermal  expansion  are  given  in  Tables  8 through  12.  These 
values  were  derived  from  limited  test  data  reported  in  Reference  4.  All 
shear  foundation  factors  are  assumed  to  be  zero. 

37.  The  mathematical  formulation  of  the  creep  mechanism  used  in 
the  computer  program  was  proposed  originally  by  McHenry."*  It  may  be  ex- 
pressed as  follows: 

ec(t)  « o Z a.(T){l  - e‘mi(t'T)} 
i-1  1 

where  ££(t)  = Creep  strain 

a =*  Applied  stress 
t = Time  after  placement 
T - Age  at  loading 
mi,ai^  “ Creep  constants 

38.  McHenry  suggested  that  a sufficient  number  of  terms  be  in- 
cluded in  this  series  to  give  satisfactory  agreement  with  available  ex- 
periment data.  In  the  present  study,  it  was  assumed  that  two  terms  would 
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suffice.  The  function  of  a^(T)  and  a^T)  used  in  the  analysis  are  given 
in  Table  13.  The  values  of  and  were  0.034  and  0.52,  respectively. 
These  creep  constants  were  derived  based  on  the  results  of  the  creep 
tests  performed  at  UC  Berkeley  as  reported  in  Reference  4. 


Analysis  Tine  Interval 

39.  In  the  analysis,  the  stress  relaxation  of  the  concrete  within 
the  finite  elements  was  carried  out  by  a time-increment  sequence.  Thus, 
it  was  necessary  to  establish  time  intervals  short  enough  to  give  the 
desired  accuracy  but  long  enough  to  retain  computational  efficiency. 

40.  A total  of  90  analysis  times  covering  26  August  1968  through 
15  February  1969  was  selected  for  the  analysis. 

Results  and  Discussion 

41.  A typical  stress  distribution  is  depicted  in  Figure  9,  which 
shows  the  variation  of  horizontal  stresses  developed  in  three  lifts. 

The  figure  indicates  that  between  days  2 and  5,  the  horizontal  compres- 
sive stresses  in  Lift  1 were  gradually  increased  because  of  the  in- 
crease in  concrete  temperatures.  The  maximum  compressive  stress  was 
found  near  the  center  of  the  lift.  When  a cold  new  concrete  lift  was 
placed  upon  the  relatively  warm  old  lift,  a sudden  development  of  tensile 
stresses  resulted  at  the  lift  interface.  These  stresses  were  gradually 
reduced  because  of  creep  and  temperature  rise  due  to  hydration  of  cement. 
At  Day  22,  the  tensile  stress  developed  at  the  surface  of  Lift  3 was  due 
to  sudden  ambient  temperature  drop  at  16  September  1968  (see  Figures  4 
and  5). 

42.  The  comparisons  between  calculated  and  measured  strains  are 
made  at  Element  235,  where  strain  meters  S052  and  S053  are  located. 
Figures  10  and  11  show  the  comparison  between  calculated  and  measured 
X-strains  and  Y-strains,  respectively.  The  calculated  stresses  at  Ele- 
ment 240  were  compared  with  the  stress  meter  Group  AG  G05  readings 
(Figures  12  and  13) . Since  the  zero  dates  for  stress  and  strain  meters 
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Figure  9.  Distribution  of  Horizontal  Stress  on  Vertical  Section 


Figure  9 (Continued) 


Figure  9 (Continued) 


Figure  9 (Continued) 
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Figure  9 (Continued) 


HORIZONTAL  STRESS, PSI 


Figure  10.  Comparison  Between  Calculated  and  Measured  X-Stralns 
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were  several  days  after  lift  placement,  these  data  were  adjusted  to 
match  the  calculated  data  at  zero  dates  In  order  to  obtain  a more  mean- 
ingful comparison. 

43.  As  can  be  seen  from  these  figures,  the  general  trends  of  the 
stress/straln  history  were  successfully  predicted  although  some  discrep- 
ancies between  calculated  and  measured  values  were  indicated.  These  dis- 
crepancies are,  however,  to  be  expected  because: 

a.  The  material  data  provided  to  WES  by  NPW  were  those  ob- 
tained prior  to  construction  of  Dworshak  Dam.  The  actual 
concrete  mixtures  used  in  the  construction  (four  different 
mixtures  were  used  between  August  and  December  of  1968) 
differed  from  those  for  which  material  data  were  generated. 
Therefore,  most  of  the  input  data  for  thermal  stress 
analysis  had  to  be  assumed. 

Ij.  In  addition  to  the  inadequate  material  input  data,  the 

problem  of  thermal  stress  analysis  was  further  compounded 
with  slight  inaccuracy  in  nodal  temperature  input  obtained 
from  the  temperature-calculation  program  as  discussed  in 
paragraph  33. 

£.  Since  Carlson  stress  and  strain  meters  were  installed  in 
groups  and  only  the  locations  of  the  centroid  of  the  group 
were  given,  the  exact  location  of  each  stress  or  strain 
meter  was  not  known.  In  a post-cooled  concrete,  a large 
thermal  gradient  could  have  a significant  effect  on  the 
stress/strain  gradient;  and  therefore,  the  comparison  of 
calculated  and  measured  values  at  exact  locations  is  im- 
portant. 

(i.  The  assumptions  and  simplifications  inherent  in  the  two- 
dimensional  analysis  may  result  in  calculated  stresses 
and  strains  different  from  those  measured  in  the  three- 
dimensional  structures. 

The  analysis  time  intervals  may  not  be  short  enough  to 
give  the  desired  accuracy. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


I 
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44.  Based  on  the  results  of  this  investigation,  it  can  be  con- 
cluded that  the  temperature  and  thermal  stress  calculation  programs  cur- 
rently being  used  by  the  Corps  are  acceptable  provided  that  complete 
thermal  and  material  properties  and  environment  data  are  available. 

This  conclusion  was  drawn  based  upon  the  fact  that  the  temperature  and 
thermal  stress  analysis  computer  programs  did  successfully  predict  the 
trends  of  temperature  and  stress/strain  history,  respectively,  although 
some  discrepancies  between  calculated  and  measured  values  were  noted. 

45.  Recommendations  for  further  improvement  of  temperature  and 
thermal  stress/strain  calculation  programs  are  as  follows. 

a.  Further  progress  in  temperature  and  thermal  stress  calcu- 
lations for  mass  concrete  structures  is  closely  linked  to 
the  development  of  quantitative  information  on  thermal, 
creep,  and  mechanical  properties  of  concrete.  The  exist- 
ing material  data  in  the  literature  need  to  be  analyzed 
and  data  storage  and  retrieval  systems  should  be  developed. 
The  information  on  early  age  (e.g.,  zero  through  three 
days)  thermal  and  mechanical  properties  is  lacking  at  the 
current  state  of  the  art.  Systematic  tests  to  develop 
this  information  are  recommended. 

b^.  Many  advances  have  been  made  in  the  finite  element  tech- 
niques since  the  original  development  of  the  temperature 
and  thermal  stress  calculation  programs.  It  is  recommended 
that  these  programs  be  updated  to  incorporate  more  effi- 
cient and  accurate  methodology,  e.g.,  isoparametric  ele- 
ments , on-board  band  width  minimizer,  non-linear  material 
properties,  etc. 

£.  The  existing  two-dimensional  programs  should  be  expanded 
to  include  thin-shell  elements  for  simulating  rigid  steel 
forms.  Axisymmetric  analysis  capability  should  also  be 
added.  Ultimately,  three-dimensional  finite  element  pro- 
grams should  be  developed. 

<1.  The  mathematical  formulation  of  the  creep  mechanism  used 
in  the  thermal  stress  analysis  program  was  based  on  the 
McHenry  equation. 5 Unfortunately,  the  derivation  of  the 
creep  constants  for  the  McHenry  equation  involves  very 
time-consuming  numerical  analysis  procedures.  A simpli- 
fied mathematical  expression  suitable  for  use  with  the 
finite  element  technique  needs  to  be  developed. 
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£.  In  the  WES  version  of  the  thermal  stress  program,  the  stress- 
free  temperature  for  an  element  is  defined  as  the  tempera- 
ture at  8 hr  after  placement.  A rational  approach  for 
determining  stress-free  temperatures  needs  to  be  developed. 

{_.  In  the  existing  programs,  all  thermal,  creep,  and  mechani- 
cal properties  are  assumed  to  be  temperature-independent. 

In  the  case  of  high  cement-content  concrete  placed  in  a 
large  mass,  concrete  temperatures  as  high  as  120  F can  be 
expected,  and  in  areas  of  cooling  pipes,  temperatures  may 
be  as  low  as  45°F.  These  extreme  temperatures  have  signi- 
ficant effects  on  the  rate  of  heat  generation  and  other 
mechanical  properties;  therefore,  temperature-dependent 
properties  need  to  be  considered  in  the  programs. 

£.  In  the  current  temperature  program,  any  changes  in  insu- 
lating properties  and  cooling  pipes  are  to  be  handled  in 
a new  lift.  Therefore,  a dummy  lift  must  be  introduced 
whenever  the  insulating  properties  or  cooling  water  tem- 
perature varied.  Simplified  methods  of  modeling  these 
changes  within  a lift  should  be  developed. 

h.  In  the  creep  analysis,  the  stress  relaxation  and  creep 
determinations  are  carried  out  by  a time  increment  se- 
quence. It  is  necessary  to  establish  a guidance  for  se- 
lecting optimum  time  intervals  which  are  short  enough  to 
give  desired  accuracy  but  long  enough  to  retain  computa- 
tional efficiency. 

i.  Solar  radiation  may  have  a significant  effect  on  concrete 
temperatures  in  certain  locations  and  times.  This  effect 
needs  to  be  incorporated  in  the  temperature  calculation 
program. 
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Table  1 

Summary  of  Material  Types 


Lift 

Date  of 
Placement 

Cement 
Content 
(sks/cu  yd) 

Fly  Ash 
Content 

(%  of  cement  wt) 

Designated 

Material 

Type 

Foundation 

26  Aug  68 

— 

— 

1 

1 

27  Aug 

3.25 

30 

2 

2 

31  Aug 

3.25 

30 

2 

3 

4 Sep 

3.25 

30 

2 

4 

17  Sep 

3.25 

35 

3 

5 

23  Sep 

3.25 

35 

3 

6 

3 Oct 

3.25 

35 

3 

7 

15  Oct 

3.25 

35 

3 

8 

25  Oct 

3.50 

30 

4 

9 

1 Nov 

3.50 

30 

4 

10 

10  Nov 

3.50 

30 

4 

11 

16  Nov 

3.50 

30 

4 

12 

23  Nov 

3.25 

30 

2 

13 

5 Dec 

3.00 

30 

5 

14 

17  Dec 

3.00 

30 

5 

Table  2 

Thermal  Properties 


Table  5 

Thermal  Conductivity  of  Insulation  Material 


Lift 

Top  Eleva- 
tion (ft) 

Conductivity 
(Btu/hr-ln. -°F 

Foundation 

940 

— 

1 

945 

— 

2 

950 

— 

3 

955 

— 

4 

960 

0.003 

5 

965 

0.003 

6 

970 

0.003 

7 

975 

0.003 

8 

980 

0.002 

9 

985 

0.002 

10 

990 

0.002 

11 

995 

0.002 

12 

1000 

0.002 

13 

1005 

0.002 

14 

1010 

0.002 

Table  6 

Average  Wind  Velocity  and  Its  Correspondli 
Convection  Transfer  Coefficient 


Lift 

Average 

Wind  Velocity 
(fps) 

Convection  Transfer 
Coefficient 
(Btu/hr-ln,2-°F) 

Foundation 

4.4 

0.013 

1 

4.4 

0.013 

2 

4.4 

0.013 

3 

2.9 

0.011 

4 

1.5 

0.009 

5 

1.5 

0.009 

6 

2.9 

0.011 

7 

4.4 

0.013 

8 

2.9 

0.011 

9 

5.9 

0.015 

10 

4.4 

0.013 

11 

5.9 

0.015 

12 

5.9 

0.015 

13 

2.9 

0.011 

14 

4.4 

0.013 

Table  7 

Heat  Transfer  Coefficient 


Lift 

Heat  Transfer  Coefficient 
(Btu/hr-in.2-°F) 

Foundat ion 

0.0133 

1 

0.0133 

2 

0.0133 

3 

0.0111 

4 

0.0023 

5 

0.0023 

6 

0.0024 

7 

0.0024 

8 

0.0017 

9 

0.0018 

10 

0.0017 

11 

0.0018 

12 

0.0018 

13 

0.0017 

14 

0.0017 

Table  8 


Mechanical  Input  Properties  - Material  Type  1 


Time 

(day) 

0 

200 


Modulus 

of 

Elasticity 

(Pal)_ 

2.8  x 106 
2.8  x 106 


Poisson ' s 
Ratio 

0.20 

0.20 


Coefficient  of 
Thermal  Expansion 
(in./°P) 

5.5  x 10-6 

5.5  x 10-6 


Table  9 

Mechanical  Input  Properties  - Material  Type  2 


Time 

(day) 

Modulus 

of 

Elasticity 
(106  x psi) 

Poisson's 

Ratio 

Coefficient  of 
Thermal  Expansion 
(10“6  x in. /°F) 

0 

0 

0 

5.48 

0.5 

0.3 

0.17 

5.48 

1.0 

0.6 

0.17 

5.48 

1.5 

0.9 

0.17 

5.48 

2.0 

1.10 

0.17 

5.48 

2.5 

1.22 

0.17 

5.48 

3.0 

1.35 

0.17 

5.48 

4.0 

1.53 

0.17 

5.48 

5.0 

1.70 

0.17 

5.48 

6.0 

1.82 

0.17 

5.48 

7.0 

1.90 

0.17 

5.48 

8.0 

2.00 

0.17 

5.48 

12.0 

2.30 

0.17 

5.48 

16.0 

2.49 

0.17 

5.48 

20.0 

2.62 

0.17 

5.48 

28.0 

2.85 

0.17 

5.48 

36.0 

3.02 

0.18 

5.48 

44.0 

3.15 

0.19 

5.48 

52.0 

3.26 

0.20 

5.48 

60.0 

3.35 

0.20 

5.48 

70.0 

3.45 

0.21 

5.48 

80.0 

3.53 

0.21 

5.48 

90.0 

3.62 

0.22 

5.48 

100.0 

3.70 

0.22 

5.48 

120.0 

3.80 

0.22 

5.48 

140.0 

3.90 

0.22 

5.48 

160.0 

4.00 

0.22 

5.48 

180.0 

4.09 

0.22 

5.48 

200.0 

4.15 

0.22 

5.48 

I 


Table  10 


Mechanical  Properties  Input  - Material 


Time 

Modulus 

of 

Elasticity 
(10^  x psi) 

Poisson ' s 

Ratio 

Coefficient  of 
Thermal  Expansion 
(106  x ln./°F) 

0 

0 

0 

5.48 

0.5 

0.25 

0.17 

5.48 

1.0 

0.55 

0.17 

5.48 

1.5 

0.85 

0.17 

5.48 

2.0 

1.05 

0.17 

5.48 

2.5 

1.17 

0.17 

5.48 

3.0 

1.30 

0.17 

5.48 

4.0 

1.48 

0.17 

5.48 

5.0 

1.65 

0.17 

5.48 

6.0 

1.77 

0.17 

5.48 

7.0 

1.85 

0.17 

5.48 

8.0 

1.95 

0.17 

5.48 

12.0 

2.25 

0.17 

5.48 

16.0 

2.44 

0.17 

5.48 

20.0 

2.57 

0.17 

5.48 

28.0 

2.80 

0.17 

5.48 

36.0 

2.97 

0.18 

5.48 

44.0 

3.10 

0.19 

5.48 

52.0 

3.21 

0.20 

5.48 

60.0 

3.30 

0.20 

5.48 

70.0 

3.40 

0.21 

5.48 

80.0 

3.48 

0.21 

5.48 

90.0 

3.57 

0.22 

5.48 

100.0 

3.65 

0.22 

5.48 

120.0 

3.75 

0.22 

5.48 

140.0 

3.85 

0.22 

5.48 

160.0 

3.95 

0.22 

5.48 

180.0 

4.40 

0.22 

5.48 

200.0 

4.10 

0.22 

5.48 

Table  11 

Mechanical  Properties  Input  - Material 


Modulus 

of 


Time 

(day) 

Elasticity 
(106  x psi) 

Poisson ' s 

Ratio 

Thermal 
(106  x 

0 

0 

0 

5.50 

0.5 

0.4 

0.17 

5.50 

1.0 

0.8 

0.17 

5.50 

1.5 

1.1 

0.17 

5.50 

2.0 

1.3 

0.17 

5.50 

2.5 

1.45 

0.17 

5.50 

3.0 

1.58 

0.17 

5.50 

4.0 

1.75 

0.17 

5.50 

5.0 

1.90 

0.17 

5.50 

6.0 

2.03 

0.17 

5.50 

7.0 

2.14 

0.17 

5.50 

8.0 

2.22 

0.17 

5.50 

12.0 

2.50 

0.17 

5.50 

16.0 

2.70 

0.17 

5.50 

20.0 

2.83 

0.17 

5.50 

28.0 

3.06 

0.17 

5.50 

36.0 

3.22 

0.18 

5.50 

44.0 

3.35 

0.19 

5.50 

52.0 

3.48 

0.20 

5.50 

60.0 

3.58 

0.20 

5.50 

70.0 

3.66 

0.21 

5.50 

80.0 

3.75 

0.21 

5.50 

90.0 

3.82 

0.22 

5.50 

100.0 

3.90 

0.22 

5.50 

120.0 

4.00 

0.22 

5.50 

140.0 

4.10 

0.22 

5.50 

160.0 

4.20 

0.22 

5.50 

180.0 

4.30 

0.22 

5.50 

200.0 

4.35 

0.22 

5.50 

Coefficient  of 
hermal  Exoansic 
106  x ln./°F) 


Table  12 

Mechanical  Propertlea  Input  - Material  Type  5 


Tine 

(day) 

Modulus 

of 

Elasticity 
(106  x psi) 

Poisson's 

Ratio 

Coefficient  of 
Thermal  Expansion 
(10-6  x in. /°F) 

0 

0 

0 

5.40 

0.5 

0.2 

0.17 

5.40 

1.0 

0.4 

0.17 

5.40 

1.5 

0.7 

0.17 

5.40 

2.0 

0.9 

0.17 

5.40 

2.5 

1.03 

0.17 

5.40 

3.0 

1.16 

0.17 

5.40 

4.0 

1.35 

0.17 

5.40 

5.0 

1.50 

0.17 

5.40 

6.0 

1.62 

0.17 

5.40 

7.0 

1.72 

0.17 

5.40 

8.0 

1.80 

0.17 

5.40 

12.0 

2.10 

0.17 

5.40 

16.0 

2.30 

0.17 

5.40 

20.0 

2.42 

0.17 

5.40 

28.0 

2.68 

0.17 

5.40 

36.0 

2.81 

0.18 

5.40 

44.0 

2.95 

0.19 

5.40 

52.0 

3.07 

0.20 

5.40 

60.0 

3.16 

0.20 

5.40 

70.0 

3.27 

0.21 

5.40 

80.0 

3.33 

0.21 

5.40 

90.0 

3.42 

0.22 

5.40 

100.0 

3.50 

0.22 

5.40 

120.0 

3.60 

0.22 

5.40 

140.0 

3.70 

0.22 

5.40 

160.0 

3.80 

0.22 

5.40 

180.0 

3.90 

0.22 

5.40 

200.0 

3.95 

0.22 

5.40 

Table  13 

Constants  for  Creep  Equation 


Time 

al 

a2 

0 

0.960 

X 

10“6 

0.463 

X 

IO"6 

1 

.920 

X 

10-6 

.475 

X 

io"6 

4 

.545 

X 

10"  6 

.565 

X 

IO'6 

6 

.405 

X 

10“6 

.500 

X 

io"6 

8 

.353 

X 

10"  6 

.442 

X 

IO"6 

10 

.318 

X 

10"6 

.395 

X 

io-6 

16 

.250 

X 

IO-6 

.289 

X 

io-6 

20 

.219 

X 

10-6 

.248 

X 

io-6 

24 

.191 

X 

io-6 

.225 

X 

io-6 

28 

.170 

X 

10-6 

.215 

X 

IO"6 

32 

.147 

X 

10"6 

.208 

X 

10-6 

44 

.105 

X 

io-6 

.191 

X 

io-6 

60 

.690 

X 

IO"7 

.170 

X 

io-6 

76 

.440 

X 

IO"7 

.150 

X 

IO'6 

100 

.220 

X 

10“7 

.122 

X 

10-6 

200 

0.200 

X 

10- 7 

0.100 

X 

IO"6 

o o o.o  n:n  oomo  oo 


_Q*tli?78.  1**4M 

I 
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COMMON/NPB  I n?nt  ap"eu 
WRITEA98, 7 0 0 1 
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400  F0BMaTUX/1AA'4> 

CALL  TRLOT 
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MBITS  10***090) 

• 05*  FORMAT  < * f PROCESS INQ  COMPLETE' » 


• 00000  AO 

000000*0 

oooooos* 

ooooooao 

•0000070 

•OOOOOfO 

•0000019 

•00000*0 

OOOOOIQO 

•ooeoito 

00000120 

•00001*0 

OOOOOIAO 

00000190 

00000140 

00000120 

800001*0 


END 

SUBROUTINE  deck 
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00000208 
80000210 
00000220 
00000222 
00000230 
90800240 
00000290 
00000210 
00000270 
00000280 
000002*0 
000003*0 
00000320 
00000350 
000003*0 
00000370 
00000300 
_f  09003*0 
000004*0 
00000410 
00900420 
000004*0 
000004*0 
•09004*0 
00000470 
000004*0 
00000490 
000005*0 
00000510 
00*00920 
00000530 
>00005*0 
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weHE»08,T99»HED 

W(rTTE»p9.?99>HED 
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1 NUNfiT4.NN00E.TlNE. tSPAN.tCARCS, JPUNCM»Soue«J>* 

2 «OtC(|>.SO|.C(3)iNPLOT.NCUAKN 

ETTg-n 


wsiTg  »8,2oii> 
58  WBtTfi  »6.2000> 
1 NUNET9 

56T 


HhDi NUPNP • KUpEU*  NUNNAT. AUPOC . NUNETji NUNfiTg# 


THJBAf«T 

OUTPUT  C0NTR0U«N0CARDiNEuCAPD;NT*PElf  -If  EACH  OP  ANY 

C ABE  » 3.  THEN  CARDS  A«e  PUACM6D  FOR  NODES  AND  ELEMENTS,  AND 

x -fuwrhFroREATEir -respectively,  *e  no  action  taken, 

RERDi NF. 2400 > NOCARD, NELCAROi NT APEli; I SCALE 
2400  raBNxT|4|») 

WITTO7.240BlNPCARP»TJELCA*D.NUPEiiJIfCALE 

....MM.**. 

C Reap  or  oenerate  nopal  point  information 
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1 ML  IN 

WBtTE17.imiN,N>0HN,X<m»VtN>t  TT ( NJ  , NUNfiVJMjRp VT,HQD# 
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S6AL8  ,'EQ  .0)00  TO  5000 


10000590 


000009*0 

00000500 

00000590 

00000410 

00000410 

00000420 

00000440 

00000440 

00000490 

00000490 

80000670 

00000(80 

000006*0 

00000700 

000007EO 

00000720 

80000710 

00000740 

00000750 

40000740 

00000770 

00000740 

•00007*0 

00000791 

000007*2 

000007*3 

oooooeoe 

eooooEfco 

00000*20 

00000640 

00000640 

00000*90 

00000640 

00000870 

00000600 

oooooa*o 

ooooo’eo 

00000940 

80000*20 

00000940 

00000*40 

00000950 

88888S58 

00000900 

0000099Q 

8OOO1090 

00001040 

00001042 


A3 


»gm  oi  n»«n-7n  i*.*** 


x»K>*nN)*i2. 
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rp-pi0Miv.i»1|i88  to  *i  - 

IF$NQM|Vt988*NBQ9Ni  AND,N«M|Y,SC,H|C  TO  «3 

gfl-Jfl  .11 

T 2 *UW<S«0HNJ «X<NJ 

TON«N«0«N|tV(Ni 

•i  ir(N.L»63.eo;7o  

» JiiJI  ijTfiil)  N 

Zf_  - i°mttly^  lUlifl  MJL 


TT 


DUl!(XiN>*g»L*l!)/OIFF 

DgMTTtN>«TT<U»l>)/n!irf 
' 80  TO  "4 


lrl"UNcy*RQ.2> 

JfiUiliLHlUfiT 


TTiU*Tl<U,l|*DF 
GB  TQ  00 

83  LL»2  _ 

•4  XiW)fiU?J»*<XONTuL»-XUN<U*i>> 
T*U>8Tax***<TUN<LU)-TWN<u;n>> 
TTtW?4>,fl 


390 


U|UU 

GQ  TQ  00 


L1«NU  

T*F*S*,l.IT0'x0}  WXll,0 
Bfo^o}  r*c 

2098>  MflB.l 


ir asv .lS»o«o}  r*crm,o 

wii«".nu‘~~  


*00  N»N*< 

Ni.N.MQd 

N2«Nt»  HOn . 

JFi<l|fpTfO>*AM»<N8,GT7&Ti  8C  TO  «fl9 
U|)T|  IO»S130> 

CatLTMiT 


493  X >N7t X<Nt I *F  ACXM  X < N1 ) -X < N2 ) ) 
T{N)*»tNlJ*FACY«<Y<Ni)"T<N2)) 

IlllLiOJ , 

IfJN.U’.HUH)  80  TO  400 


l*W 

80  wS}  T|  to , 2002 ) tK;x<K>,YjlK|#TTcK> , K»Ui*tT 
l*L*l 

_ -- 

jFjMoOtQT.0)  00  TO  890 
IP  (M"t>  to* 80*79 

9Q  IP  180.100.80 

100  CONTINUE 

IFJNTIPElt,EQ.0)8O  TO  3qOB 


80001011 


80801090 

80001000 

80001090 

000010*0 


000010*0 
_ 80001111. 
00001110 
00001120 
iOJOllll 
000011*0 
00001190 


000011(0 

00001148 

00001180 


00001210 


00001220 

80001290 

00001200 

00901290 

000012*0 

00001270 

JOOO12O0 

’ 000012*0 
00001310 

80001310 

00001320 
00801390 
. 00001300 
80001390 
■00013*0 

S5001320 

00001380 
000013*0 
_ 00001400 

000014*0 
00001420 

00001490 

000014*0 

00001490 

000014*0 

00001470 

000014*0 

*00014*0 

00001500 

00001510 


J279T  01  09»11”78  16,496 


W»tTSUl»t»*K«N)*r(N>»N«l,NI.M»iP| 


3000  IP«N0B4R0»eQ.0}00  TO  2«10 

ptificsr  b #xm>TYTB  »;tt< 


2*00  rOBH*Tu5l5x;3FX0.3> 


X 8EaB  and  PRINT  of  ELEMENT  PROPERTIES 



3610  Mtinil.itHl 

N«f 


•0001*10 
90001919 
•0001990 
00001990 
90001900 
00001920 
101910 


103  PRADINF  >10021  H» ( I*(M<  J ) * I»1 •*  > *PLT(i(H)  *HOB*miN  iTi  MPlTH*  PtTlNC 

WnTB»7.xf02|N.(JX(N(nlI*ii9),PlTH1Hj,N0u;NtlPlI*H9tTN, 

1 BfeT JNC 

IF  (MOO  .N«.  0)  BRIJ|(Of 2099)  NODfNUHIT|NeLTH«PtTlNC 
__lfi6  _NSM  ^ J|_g-  X07iiOT[l05 
109  IR{N.l»*JX|N-l,ll«l 

IJLJNf  2 ) !I8|N;1 1 21  *1  ... 

!R{N*B|»!>UN»l«3l*l 

IR{N.jl*IRlN-l,<»*i 

_ _!«*' •l*isn5i,?j 


GO  TO  107 
106  IFJN.BJ.NUIMIT}  90  TO  108 
N*N*1 
NlaN.MOO 

n2;N1»hqb  _ _ „ _ 

TTJTn ftU T .0 > . *KlO»Y N2VGT70Tr  Ge  TO  110? 

mMt6  16.2130) 

CMjLjXlT 

~U0T-Dirn97  jj|V,R 

210?  !IJN.I>»t»IX(Nli|)-JX<N2»I) 

— iorUlTe’ii  *’*<  Txor,T57rii75  jtpltm  jn  > 

IFJmoD.Eq.O)  go  To  6107 

}FtM.N6.W>  on  Te  106 

tloMS 

RRiNtINIT*MPl,TH 

_ 09  9107  mi  .KKIhPlTM  


DO  4107  J»1»MPUTM 

mTl*i  V »pl  tm  t i.r 


tist  _ 

910?  Pu.THa*lJ»9LT«UTl>*PLTINC 

108  IF  (NUMEu'N)  10*il0*il03 
10?  fflTlNPE  

US  Pt.TMU>«PUTMtI)#a4, 

! DO  119  N«l*190 
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0000199b 
OOOOUOO 

90001610 

00001610 

simiii 

00001690 

00001690 

00001670 

00001600 

00001690 

00001J90 

00001*10 

00001730 

000017JJ 

00001760 

00001790 
00001790 
00001720 
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000017*0 
00001000 
00001810 
00001820 
00001820 
00001860 
00001890 
00001860 
00901870 
00001890 
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00001*90 
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000019*0 
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00001*78 
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IP{NUM6Ti  ,E0,  01  GO  TO  121 

ISlSi. 

rbac 


*0002010 

00*02020 

*0002000 

00002000 

•0002110 

•0192110 

00002120 

•0002110 


•021jD 


8 40  < ff » 1 0 0 • f t (IT ( 1 1 J*  i j * J *1  # 2 ) , ! *1 1 HU^E  ! 1 ) 
«2tTE47.l008H{BT(I,J.i)iJ»l*n,I»i;NOM6TlJ 

• 2)H»1«N0KEU) 


MiiltAiil9121_IET. 

0®  112®  i*1.MUmFTi 
UZO  ET?I.t»l)4ET<l.i;i)»24, 

_m_H.SNimfl2_.lfli  01  Gfl_rp  JJ2.  _ 

«B#0<Nr.i|08»  UET(1,J,2),J«i,2>,!»i,NUN|TJ) 
NRlTfilT, liPOK tBT ( I , J*Z ) , J»1 ,2 ) i I»1*NQHET2> 

! • 1ET,<(ET(j,J,2)(j.i,2),l.},NBpET3) 


“RTTEIO. 20121 
D0  1121  I *1 1 MUMET2 
1121_ET{I,l,2)»ET(I*li2)*24, 

122  HtNUHfTS  .EO.  01  GO  TO  129 
IET»J 

»B4P*HF,l*081  t(6T(I.J.3)>J«1.2),I<il,WUHETl) 

URltEiT,i009MfET<  I.  J.3T7j»l.i)»  l«iiN0MET3} 
«B|TE»*.20121  IET,<<ET<i#J.3),J.i,2|,:»l,NBHBT0> 
DO  128  I«1#NUMET3 
1 23Ft  t H 1»3) » 8f(!  i i I S ) • 2 4 , 

C FOB  6ACH  INTERVAL  of  T 1dE»*SEVeR At  TIME  S!|PS 
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•00021*0 

-ttlttffl 
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•0002210 

•0*02220 

•0002221 

00002210 

000022*0 

000022*0 

000022X0 

000022*0 

100022*0 

00002300 

000023*0 

•0002320 

100023*0 

•0002340 

00002310 

00002340 

00002370 

009023*0 

000023*0 

00002*10 

_|00024*0 

00002420 


00 

00 


243J 


125  CONTJUUE 
NcgHfO 


244 
000024»0 
000024*0 
•0002470 
00002480 
000024*0 
•0002500 
00002510 
•0002510 
•00025*0 
00002540 
00002510 


A6 


ooor 


327*T  01 OjtfU'78 

N«SHM 


190 


133 


Da  130  Njl.NUMNP 

Htl>*l«0 

77HClBDIr»D;fl)  "ID  TTJ  180 
KC0RO  «Ie«KD«/9 

jg  ■«CArD».  KCjjjPO  

If tKZ  ,Ne,  01  KCARDpKCARD*! 

Kc  ■ 9 

Dti3*  UlrKCARP 

jf|»*  *eqt  kcapdi  .and.  «kz  ,ae,  on  kcmz 

R*80fNF. 1*391  tl0(N)|TA(N|«ktl«RC) 

“11^1 7 f J"iN  JATR  « Ji  I iMHi  *C  ) 

* ■ :«iJ)  ' 

TIJO  * T«  IJ1  

CONTINUE 

PffBHif  <a<i9;rlO,9)) 


c 

c 

c 


13*  iPjNNOPE  »EQ»Oi  00  To  140  

«S*D<NF,l*10)  (N0DE( { ) , {al.ANCDE ) 
wmTE»7,l91CMN0DEU>#  m.NKOCE) 

_ RE*D<NF.i0i0  » (|.IFT(l)i  I ■!. KNCDE  ) 

«« tT E * 7 , 19 1 0 M U F 7 ( l ) , l ■! . NA OC E > 
mbits ii.iaao » 4node(I).i*s,knodei 

WHITE  46. 1*01 » tLlFT(I), IPlihNCDE) 

ir  up^ncn  .RI.  onrtiTe  ,6,2035)  inuncm  ~ 

IPUNCW  « 1 PUNCH  • 24 
ReaF  AND  PRINT  OF  LAVER  PROPERTIES 
l|?AN  ■ !0R*N*1 

reading.  90«.  end*?97itjti.e  • 

Wf»ITE»7,9*9)TITl.E 

RE*DfNF,i*0«)  NUMN.NUHfc,NUMCB,NtMCPilVpr,  INTER, CT#PtTINE, PLACET 

■R«rTr(77i0a4)NOH1l,THmr,llUHCt1NUPcF,BDT,.INTIRlBT,PLTlNE,  PLACET 

IF  (NUNN.aT.O)  00  TO  191 
IF“(NRLOT.ffe.OrSTOP  “ - 

RrTURN 
CONTjNUE 

MB  JTe  4 6*  2*32 1 I SPAN , T I T(,E 

WB1TeI*.2*05»  NUHN,NUHE,NUMCB,NI.MCPINBt*  InTER, ST jPLT IMS, PtAC«T 
PLTlN8fRLTlR5«24, 
rFtTSRAN.EO.l’NEtP"! 

ELIMINATE  all  INSULATING  ELEMENTS  AND  CO0L1N0  BIPES  HJTH  REMOVAL 
TWE  LESS  THAN  PLACEMENT  TJPE  AND  READ  ADDtTJOlAL  ONES 

KNfO 

DO  159  N«1#NC9H 
IFfTMCIN) .LE. PLTIME)  GO  TO  159 


in 


10*0291* 

•0*02929 

000023*0 

•00023*0 

•0002*90 

•0002*10 

•0*02*20 

90002*00 

00902**0 

00002*90 

00002**0 

00002670 

00002*00 

00002*99 

90002790 

00002710 

00002720 

»;n 

00002790 

00002760 

00002770 

00002700 

00002790 

00002000 

00002060 

00002020 

00002000 

000020*0 

00002*90 

90002060 

00002070 

000020*0 

000026*0 

00002900 

00002910 

•0002920 

00002930 

000029*0 

000029*0 

00002960 

00002970 

900029*0 

•000299* 

00002995 

00003090 

00003010 

GO0030J0 

•0003030 

000030*0 

00003050 

00003060 


A7 


JHZU.JL.  A IslizZi 

KHRKK*) 

JC«N> 


it: 


JGfKK’fJCt 
HCtKK)fHG(N) 

_TmC«kK>»THC(N) 

— ct  jinri  

I6TC(N 

JB-mtth 


fcCfjL*’ 

NCBH#KH*nUmC0 
IFtNuHgB1jOJj)JL^0_TO_ 
NffCBL 


1*9. 


6*6  *i*D<NF,lfOSI  tC4N)(JC<N)»HC(M.TFC<N),|BT|(Nl«IMN 

IFtUGN  *10.  0}  90  TO  IS* 

JJ«<JCIN).IC{N})»1 

JClNJ»lCiMJ*Jk  - .. 

08  1ST 

1C{N*J  >*IClN)*l 

JCtNU)»jClN>«l 


HC{N*J)«hC»NI 

THCtN*l>1THCtN) 

I£TCiSni»XEICiiU ...  . 

VS?  N*tl*i 

isa  n«n«6 

lFtN  *LE.  NCBH)  Bo  TO  166  _ 

WBfTEli.iyiS)  l jC(N),kC<  N ) ,THC<  N > * I *?C < N ) * N i Mi 


I?? 

160 


DO  15?  N«NCBl,»NCBH 

Td£<N>_fTflBJM»»Z44  _ 

K|MO 

DO  169  N,l.NBpH 
IF(THI»|N»>LE>PLT1WE>  ttO-TQ  165 
KK*KK4l 
!PtKK)»IP(N) 

HPJKK»*HP{N> 

?e{KK>Pi»ii») 


• NCBM I 


1*5 


Tgg<KK>*TMp<N) 


_ SriHgjE 

NCeL?KK*l 
NCPHfKK*NUMCP 
If  JN^MCPaEQ,0  ) 00  TO  170 
N»NCPt 

166  RiFDlUf .19061  IP<N>,HP<N),TP<N),TPP<N},V.ISC 
WBtTEt?.H06HP(W)>HP<N><Tp(N»«TMP<ll>»JMNe 

N«W*1 

If  {N  .GT.  NCPH)  QO  TO  168 
IF  UNO  •GO'  9>  00  ?0  1*6 

167  lP{Ni»jP,S-l»*lNC 

HP$Ni»MP(N-l> 

TUN  IflRi  l|=l.|  .... 

TN?(N>*Tm!!|%1) 

N*N*V 


•0*911*1 

•00039*9 

*00*3160 

•006)110 

••••91I0 

iiiiilii 

•0003160 

I0009WO 

•00031*0 

10103101 

00003911 

00003260 

00003230 

00903280 

000032*0 

•00032*0 

00003260 

00003270 

00003200 

000032*0 

00003310 

*0003369 

.69603320 

00003330 

000033*0 

009033*0 

*00033*0 

•0003370 

•00033*0 

000033*0 

00003400 

000034*0 

00003420 

00003430 

000034*0 

000034*0 

009034*0 

00003470 

000034*0 

000034*0 

00003590 

000035*0 

00003520 

00003530 

000035*0 

•00035*0 

909035*0 

00003570 

00003580 


_32.79T__p.l__  OliilrlB 19.06  

JPtIP(B-ll  .LT>  J>  90  TO  1*7 

1P6N  ALE.  UtfHJ  10  TO  .1** 

163  NR[TE1(. 2(141  U^«n>|HP}N),TPCN),TMP<N>#N,|«NC2U;NCPH) 
DO  16?  N»NCPL*fCPH 
_cl*5_It)P<il»lTJlMN»»2«l  . 

6 CNfiCN  lNCOHS*STf|(CY  Of  tlfT  jAfCRPAflON  a*D  RgpOVAl  Of 
_c LJNeNIS  IflD^Og  900WLN6  PJP|« 

170  KM 0 

VY*DTpFLOAT(NDT) 
lFtNCB(J.eO»0|  96  TO  261 

D>  209 

**»! 

200  COST  (NUE 

GO  TO  20“ 

DO  20*  N*1 »NePy 
XMTMP*N).PLt!ME 

5?iXI'tT.I»>_KH»l 

209  CONTINUE 

208  JFJKKJfjE.P)  HRITE<6,2019) 


“C "SW^AU,  NI*  NODES  To  PUaCBNINT  T6PP|RaTURE  ARC  CONTACT 

C A*6RA«fi  TEMPERATURES 
C 

urrn  Tir*Twn»  ~ 

en>»o,o 

210  otmi.o 

DS  221  Nal.NUNf 

jFJPtTHlN) .QT, PLTIME)  «0  TO  220 
OR  2(9  1*1*4 

Ti*mN.n  _ - 

IP  fPLTN(  N)  .f:Q.  PLTIME)  b ( 1 1 ) »B  (I  I ) pPLACET 
I FJPLTM(N).LT. PLTIME)  IM  1 1 >«B< 1 1 >*T< It) 

219^0TU  Tf  BIT  71*1.1 " — 

220  CONTINUE 

DO  238  Nat, NUMN 

rPJ9fNT.69>0»irT  "90  T 0 '230 

T»N)*B(N)/0(N) 

230  CONTINUE  


CALL  LAYBR 

IP JTeHPMaA . GT .pgakt§mp ) 

TgltPHA  i .Pg  *KTeNF 

poattpsakDay 
NPgAKNDEaNOnEPEAK 
J766  G'9  TO  (4f 
C 

C TOIMaT  statementq 

?97  HRtTE»(*9T8>  * 


>Q0  TC  *700 


•0009991 

oooo9«io 

•00036(0 

00003620 

00003690 

(N9UL0tIn9  00003690 
(00086(0 
•0003670 
•0003600 
•0003600 
00003710 
•00037(0 
00003720 
•0003730 
•0003760 
00003710 

oooos'te 

00003770 

00003700 

00003790 

00003*80 

000036(0 

SURFACE  AT  00003020 
00003030 
00003*60 
00003010 
0C703663 
00003*70 
000038*0 
00003890 
90003900 
000039(0 
00003920 
00003930 
■“  00003960 
00003910 
000039(0 
00003970 
00003900 
00003990 
00004090 
000040(0 
0000*040 
00004030 
00004040 
00004C*0 
000040(0 
00004070 
000040*8 
000040*0 
000040*2 


L 


t»l  r««HAV|*iN»  er  nle'i 


f*Bh 


9*9  F*BMA*»16A4> 

IQOp  FORMAT  (•I9.F9,043I»«3F9i0,!2,J3» 

iia  

1101 

tool 

1004  FORMAT (6 | 9«3F10 , 0 ) 

1009  FORMATt2|9»2Fl9.0.2i5> 

A P 0„‘_FRfi"AT  u 3EU. 

1900  Flf»»ATt2rl0.0) 

FfRMATtlA»9> 


FORMAT  <4t», CIO. 9, 3|»,F9,I> 
nlMAT  (t»»Afl*iB> 


9090*0** 

-101040*9 

90904190 

00004110 

99904120 

90904190 

9010*1 *0 


i? .iulutra  auTMJ'/iMiln 
out*i 


|MAT114tO> 

jTLLaMOI „ 

1031  F0RN4T|'  COLUMN  HEADINGS  FOR  extra  0UTPUT»2(24T5|) 

2000  F8R«AT'('9*.  SX,16A4//25MoNUM8ER  OC  MOBIL  POINTS*-, U/ 
1 2»M  NyHggR  qf  ELEMENTS-——-  J4  /2»H  NUMRRP  CK  MAT|R1AL5— 

2 JJH  RRHiRR  OF  QI  cards——  14  / 

3 *»5H  80  OF  S*T  TEMP  CArOS—  |4 > > 

_299i  /JfiMtLLlilllL  NiPx  Mfl1^4jc1ifeljl4lt,lilT*llK*#HTSeP| 


2002  FORMAT  (i|in;3E10.6) 

2003  FOBMAT  (51H0  N I 
2904  FOBMaT  (SiSjJIOjEI*,*) 

I860  PflNlT  <19 1 /'JOHONuMBfeR  OF  NOPAL  POINTS 


L MATERIAL  PLACEMENT  TIME) 


14/ 


00004140 
000041TO 
00004190 
00004190 
00004200 
- 00804210 
00004210 
40004290 
14/90004240 
00004290 
00004240 
-9090*2?0 
00004200 
000042*0 
00904319 
90904310 


1 >9H  NUMoeP  OF  ELEMENTS—*  |4  / NUMBER  OF  CONVECTION  PC-!4/f non43(B 

2 «»H  BUMbIR  OF  COOLING  PIPE-*  14  /J5R  NgN|£R  OF  1NCREHENTI,*— 14700904388 

29H 


3 2§H  OoTput  JNTeRVAU H / 200  TtM|  INTERVAL——  F10.3/  00004340 

4 29N  BEGINNING  TIME FB, 2/JSlH  PLACEME8T  TEMPERATURE*—  00004300 

5 FJ,2>  _ _ 00004310 

~3006~niMAV|iM9  M,11X,4hC0ND,U*,4hSPhT*UX,|hcENS  ,4X,2*hTIME  MEAT  0000432Q 

- - --  _ - — - - 00004300 

TIME 10004340 

1 1X.UHTEMPERATUR6J  00004400 

2009  F0BMATIF9^2,E15,A)  00004410 

m?  FORMAT  (4H0  M J4X  1HK  J4X  1HC  14X  JHD  14X  IN®/  <]4,4E15^6>>  _ 9000**10 

1 FOB'MaT  T27H1TW0  DIMENSIONAL  PLANE  BODX  I 90004490 

2012  FBBMaTv27U0TemPe0aTURE  OF  ENVIRONMENT, 12  / 7H0  TIMe,4X,  40904440 


*A»j*  

IQENERRtlON  STOPS/( 14, 3F15V4,F30, 5)) 

2907  FpRmaT|4|UQA0I9BAT1C  ThRPERATLRE  RlJE  Of  The  NlTgRlAL/  9M0 


JoIsS-  O^M^fuofN^ut  T NG^fc  LEME  kW/  4h 
IEMOVeo  ENVIRON  aRRaT  POSITION  /( 


t,*M  «',1*X,1HV,*1N 
2!?.*ri3,«,T5)iio>> 


_JI00044fO 

TIME  T«00004**0 

■ w»,  ' i»"  r \ A i 7 , *,  i]  | o , i 3 , t ill  I 1 00004470 

2014  F8BMATX25U0DBT AILS  OF  COOLING  PIPES//AH  T , 1*8 , lHH, 14H  T|MPERA000044}0 

‘ iTDBE  *9oh  TIME  REMOVED'  ARRAV  PoSlTl0N/ll5.3Fl5,*»U0>>  000044*0 

2910  FOgHAfilM0,4X, IMMATERIAL  TT*E  *,13)  00004990 

2019  F9BHAT<39Uo**»»  ERROR  QESS* Of  _ ..WARNING  OJLX/IX; TlWNEW  (JfT  1ATA  00004510 
lis  noTsubplied  even  though  a change  uas  oicurid  in  Insulating/  oooo*?2o 

■■  " ALeULATtON  PROCEEDS!  00004590 

00004540 


2 '/X,51HELEMENTS  AND/OR  cOOLINo'PIPeI,1  PaLCULaHoN  PROCEEDS) 


2g|l  FBBMaT 


(ISHOBAD  CARD  NO.  14> 

}32  FBBMAT^  ' ISP AN  ■ ',li/X*A4> 

2035  FBBMaT  USTeHPERaTUReS  PURCHEC  AT  p^Y  • IJ,  l FOP  NIITrRT  Rj|"t 
20*0  F8BMAT(1H«,60X,5H  MOD*, I3»3* 1 *H  NLlM% 14, 3* , AH  FACX*,FlO,5, 

1*8  FACX*,F10  »*) 


00004590 

900045*0 

90904920 

00904540 

0000*5*0 


3279T  OX  B»t11«78  J.*,496._ 

3999  P b8NA*  * ' |60X , • HODi ' » ! 5, • NL JP !▼-* , IS# » HPlTp»»,l»,»  PLtJNC»'. 

^T^^hat*42no  insufficient  infcrp*tjcn  to  6eN|R|T|  *eSh> 

END 

SUIRftUTjNP  LAYER 

PaBANBTER  HANN95*0*HAXBR«9,KXI«<eo 

CvPHON  NODAyiNHHfil.NCB^NePMUPHATiNDT.JR.TEMTiTJME.NUHNPiNUME, 

1 NUN9C*PLTIME»8fHAXN> #X*NAXN)|V<PAXN|#T4H|XNiIDInABNi i7T<HA|Ml • 

Z TXIHjr.9y,FUTR{MXE).VOUIHXE>,»-EC<i4j,CP«5),S{5;3)»KX(«>.8(S#5). 


gCONP  < 10 ) 1 8PMT<10 ) (DENSIXO )iOX(35J2*10 ) iNSTOB (10 ) «ET (100|2«3 ) • 
|C<1JOI/JC<Z10>.HCI210>*T^C(210),CI.C210,4ISP#A,PUACET,NELP« 


3 

4 

"5 

6 ?SAKBAY,NODj|PEAK,PeAKT6HP,APlAKNCEiT|F5HX|> 

i'SOtCf 3 ) , NPLOT i NDaTI <2l 
COPMON/SVMARB/NUMN»  MBANDi  A (I>AXN*hAXBW)»0(N|XR) 
COMMON^NPBIN/nTaPEU 
_T)IMEN1IOn  TnoDE<10> 


10l<.JC<210».HC«2AO»»TPC<2ig,*CUl210»iISP*N.  _ . . _ _ 

lPJli0j,HP(210),TP(2IO),Tt<P<310j,IETC{319)lNE,TEHPMAXlPDAY, 
_S*KBAY,N0D|PEAK#PEAKT6HP,APlAKMCE»TfF jHX|) 

CBpHON  /OUTPUT/  NNOD|#NcOUNT,KOCE<30);Ur!t3B|iTiTt|(j*),jPUNCH, 


ToRfr  BBWCTTV 1 Tr>» At RTirr OB  COMPLETE  body 


130 


DO  130  !»1iMUmn  

~D*tTo*»d 

BU>»».0 

Ott>»B»0  

isoi3o^r*iywA"xBw 
A‘t.J>«0»0 
MNAND'O 
TSTWtl 

_OJL.ll*  Nltf  NUME 
t F {PtTN ( n) , qV.PL TIME I 00  TO  200 
MTYPB*IX<Ni5) 


Cfl|p»2C0wP fRTYPfi 


ELWffT  CONDUCTIVITY  HaTRjX 


1»0 


DB  IBB  J*i#5 

"infTT>TX(¥,rF 

Do  191  J»1  • 9 


_k*UP<IT 


j£ms± 

EUl*lUl*ULU/li_ 


YI1(TU)>1< 

vijtTfiJ*0j 


0000*800 

00004410 

00004430 

00004430 

0000444B 

0000««4B 

00004490 

000044*0 

oooo*«“o 

00004400 

00004400 

000047|| 

000047(0 

00004730 

00004730 

00004740 

00004790 

000047*0 

00004730 

00004780 

00004790 

0000*800 

000048*0 

00004830 

00004880 

0000***0 

0000*890 

000048*0 

00004870 

00004840 

000048*0 

00004900 

000049*0 

00004920 

00004910 

00004940 

00004910 

000049*0 

00004970 

000049*0 

00004990 

00005000 

000050*0 

00009030 

00005030 

00005040 

00005050 

000050*0 

00005070 

000090*0 

000050*0 

00005180 


All 


n 010 


»27*T  01 0»sllT7» 


Dt  193  M«< 


90009110 


"7»LHiKr 


■OOOglgO 


If  (I’lJ)  135*192*139 

i»»  *4«xiUI-xtl) 


A**X*-*<|} 

fMJUUL'jrilL 


BK»Yy«*ui 

C*|J-0K 

D*?AKyAJ 


XL*M**J»BH.AMiJ 

If txL>B« 9Tif0>0 1 to  TO  Iff 


IBT0P*1 
U?tTE(«. 29031  N 


__  i?6_v9k(ii1iVQi.«  n ueuhmm. 
Cl8H»*>*c* 


•oB§d/xlxh 


g|iil)lC.*»t*JX»gl_ 


90009110 
00009100 
00009190 
00009100 
00009170 
_ 900091*0 
900091*0 
00009299 

mKiH 

9000929Q 
00009210 
_ -00003270 
9000*290 
909092*0 
90009310 


:*1.2»»B 
E«l,3)*-BJ«OAJ*DX 
_ 

E»3.2)fBnf»2»AK**2 
El2»3>*-Bd»8K-AJ«AK 
Et3»l)fE(l»3> 


90009310 
30053*8 


90005330 

00005330 

90005300 

00005350 

90009300 


E}9,2)»E(2*3l 
:83.1>*‘ 


•Bj»»2»AJ»«2 


Mill** 


^f*I«** 


140  kX{2)*1 
145  K9{3|*> 


M)  t491140;i49 


DfT91  (>1>3 
II’KXIJ) 

DO  1*»  Jilf3 


00009370 

00005390 

..00(1053*0 

000054(0 

90005410 

JIUlill 

00005439 

00005440 

.90005490 

900054*0 


900054*0 

sooo949o 


J4*KXIJ) 
191  S»tl. 


UJ)*9<I!,JJ)*E(l,J>*C0fM 


i9a'c#wf ji«E‘ 


39  199  I»l,4 


900054*0 

00005590 

00005510 

00005510 

90005530 

_Afl0055ig 


DO  155  J.1,4 

155  S*Tr J)»S(lr J|-S( J*5>»S(J»5»/S(5.5) 


J,  90D  BUEHEnT  CONDUCT  I VITY  TC  COMPLETE  COWDCYlVITY  HAT*!* 


VOtlN»»VOLlNMSPMT(WTYPt)»DENHPTYPi)  90jJ? 


DO  179  L*l* 4 
I*LH|U 


00005550 

000055*0 

00005970 

000055*0 

000055*0 

_0.0005«80 


000056*0 

00005620 


A12 


3279T  OX  Q(«ll-»7B  ift.456 


Dtt>*OU)*VOL«N) 

Pg..i7»  "ill* 


_1?7 

»l 

L*l 


}!? 


IF  (HAXBW-J)  157*i50|1»( 
WHITE  16,2002)  N 
00  To  300 

!FtNB*SD-J>  16?, (65, 165 
HgjWp.J 

iMjj  i75ayr.iTo 

*•». J»*A( 

CONTINUE 


20fi  CONTINUE 

IMloy^H.* > STOP 


BOUNDARY  CONDITIONS 


TPWW.80.0T  SO  TO  220- 

DO  3(5  N»lrNCBH 

1»1C(NI 

J*dC<NI 

Xl.»S0NT<lMJ)-X<I>>»«2*(YIJ).Y(I))»62) 
H»4lC  <J<I  2& _ 

IKJ(K?Sr2T2T273To  ~~ 

A*lMO»A(l  >Kt*U 
_s?  TO  ?15„ 

Aj J#KI*W 

CL 


«yNim- 

^DNTjWE 


COOLING  PIPES 


IPtNCPN.8fl.0l  GO  TO  225 

DB  22*  Njl.NCPH 

i"*ip<n) 

Att,l)pAiniI*HP<N> 

22A  B»t)?BU’*WP<N'«TP(N> 

225  CONTINUE 
c 

C 2r  TBMPBBATURE  BOUNDARY  CONDITIONS 


00005«>0 
• 00096(10 
00009090 
00005600 
00005670 
00009600 
•0009600 
00005710 
•00097(0 
00009720 
0o0o97so 
00009700 
00009750 
•00097(0 
•00005770 
00009700 
•OOQO570O 

000050(0 

•0009010 

00909020 

00005810 

00009810 

00005850 

000098(0 

00005870 

000058(6 

00005800 

00005*00 

000050(0 

00005920 

00005910 

00009*40 

00005950 

000059(0 

00005970 

00005900 

00005908 

0000(090 

0000(0(8 

0000(020 

0000*030 

00006040 

0000(090 

0000(0*0 

0000*020 


DO  300  Nal.NUMN 

C 

0000*000 

OOOO(0y9 

DrmNT  Po.e.or  go  fo  sob 

D0  250  H»2, HBAND 

K*N*M*1 

0000(1(0 

oooodio 

•000(120 

— TfjKJ  235.235.230 

230  B»p)6B(KJ-A(K.M)»TTIN> 

0000(120 

0000*140 

J27*T  01  m»ll?78 


AlR.MtfO.f 

***  inK!«5-t»  545,2-40. 240' 
240  B»L)*BtL)-*(N,M)«TT(N) 
249  *>tl*H)»0A0 

TlfiHTlW 

ruumiMi 


300  CBNTjNUE 


lOIWT  7fW|"6T\JRE* 


ranm  wftroffammtiai \t%mwnmr 


J><L  3«0  W»1,NUHN 


•000*1*0 

>0»0«1»0 

•000*120 

•OOOtllO 

•00061*0 

loitiih 

•0006210 

•0006220 


•0006230 

•000*2*0 

>•00062*0 

0000*260 

>000062X0 


!^t*ll»H>*oo»0*0J  Mn>1)»3«0 
jrjTnBj.He.f.*)  G°  TO  221 
PMI>*BT2»PTN» 


320  COWTiNyt 

c*u  *rHaBt<t> 

sbt  o*  Fia*s  for  stress  free  Terr  outrui 

— - 

NCriG»i 

NCSUI ....  

II  THU  a fake  llFT,IE.FOR  R 1*8  R6H0V»l  IN  HIDLIfTT 


iTTNELR.QT.NUME) 00  TO  7(0 
SFT»8? 

KOT«SFT/pT  _ 

JFJDfvQt.SFTMSO  TO  702 
IFtDT)60.SFT»Q0  TO  701 

PPT«HDT 

DUTltSFT/PT 

D6l«DDU-PDT 

IFfDBUGe.0,5)G0  TO  702 
iFJUBltEO'O. )G0  To  701 
700  NC!l6NpT 

NCFl  11*9 


Go  TO  70s 

701  ngil»kdt 
ncflg*4 
QB*t8  705 

702  NCT16RDT.1 

NCPLSfS 

705  CONTINUE 


ssisim 

•000*320 

00006330 

1000*360 

0000*3*0 
00006360 
_ 00006370 
8000*3*0 
•09063*0 

100063*1 

00006400 

•00064|0 

OO0044J0 

•0006422 

00006424 

•0006426 

00006421 

00006420 

80006460 

•0006490 

•0006460 

0000*^0 

00006410 

000064*0 

00006500 

000065(0 

00006520 

00006530 

00006540 

00006510 

00006560 

00006570 

000065*g 

000065*0 

00006600 

00006610 


A1A 


7 ^ 


3279T  01  Q»iU-78  It. 454 

C C*tC  7«M»  AT  END  07  EACH  TIME  INCREMENT 


C 

C 

T~ 


determination  or  heat  uengratiok 


DB  3*9  N|t«NUME 
irtPerM<N>.8T,TIHE)  00  TO  299 
~ mTjPeA jX(Bi5| 

T*»T jM6-PLTM(N> 

lrjT*4SE.EST0P<HTYPiJ)  JO  TC  2*9 


384  D!{T»aXcul.NTrPg,-<Jx<l.*l»l.MTYPE) 

GRAD*i8XlJ,jL2:MTYPE>-0X«l,-l»2»PTYPE)J/DIfr 

OBi8R*ppyou«M» 

DB  391  l«l»4 


"5Vo  0 

39*  C 


S5Wn*0Q 


COB VEIT  ION  BOUND ART  condition 


I«tC(MI 


TO  410 


1 


*jci Ml 

mi'i 


IIETG(M) 

dq  <g«  n*i»ioo 

H-N 

xz*ET  (n#  j,»  i 1 1 "Time 

IFCXZ)  4Bo,480r401 
400  COfiTjNUE 

fOT  VirffSTf iJtllllTrETf  I J-fit,  IV 
XJZ»XZ/ DlCF 

Xgjt*ST4l  J/.2*  l 1 >»Ex<  l J-1.2*  1 1 ) 

Tjjnp_?rn  j j ,211 1 y.XXX  * JtYZ 

TeNP7T6MP*SOLCUI) 


XZ«HC«M)»CL(M>*TEMP*O.S> 

*n 


c 

c 

-C. 


UitHDU) 

409  Q»J>*81J>*XZ 

iio  cosrjNue 


lv  CAUCutATB  EFFECTIVE  LOAD  MATRIX 


“WO 


“4  91“  J a!  # Ni/NN  ' 
0«t>l8U)*B<I)*D<I>»T<l) 


lOQOAAII 

90006616 

00006426 

00006630 

•0006668 

•0006690 

•0006660 

000066^0 

•00066*0 

•00066*0 

00006780 

000067(0 

•0006720 

80006720 

00006760 

00006790 

00006760 

00006720 

000067*0 

*00067*0 

00006*80 

00006860 

•0006620 

00006630 

00006660 

000066*0 

00006860 

00006870 

000066(0 

00006690 

000067*0 

00006*10 

00006*20 

00006*30 

00006940 

00006950 

00006760 

00006970 

00006980 

000069*0 

00007C96 

00007040 

00007020 

00007C30 

00007060 


000070| 
00007070 
00007060 
000070*0 
00007180 
00007140 


327*T  01  Q0»U»78  It. 496 


2*  SOUVI  FOB  TEMPERATURES 


C*UU  SVM9BU?> 
DO  5qg  l«lfNU*N 

9oq  ott>t*«o 


TjpE  «TImE*DT 
DaI«T!(1E/8*, 



PUNCH  CAHD9  FOR  OESTART  run 


1fTJJJ,NE. JPUNCH)  go  TO  5*0 
WBtTE«97,2004j  <J,T(t),l.i,AuPN) 

LFiHfifiM  *|Q.  0>  . 

Do  sq»  N,*,scBy 
XTPC  » TrC»N|/24,0 

5 91_ WHTTE1I7,2P0Z>  1C( N> » JC«N)  • * T|)C* Jfe TC { M 

502  IFJNCPH  .M,  0>  qo  to  510 
Dp  503  NtJ » NCPH 

XTMP  * T^E«NJ/Z4l0 

503  H«*Tg»B7,8O0«)  JP(N)»HP<N)*TP<R)#XTrP 
2006  F08maT  <5U5;f*0,5>) 

to  oe’foBHAT'l  ?5?5x*  fio  ! sJriolifri  o «~ i > 

510  iFfNNOBE  .EO.  0)  GO  To  954 

_.LP{ABQP(nAyji.9.i..*Ne*  60  tcjs* 

irtNODfid)  • •*•  nunni  g6  to  ss* 

xegToi 

TNQDEM>  • 0.0 
lt«N0DE<|> 

jFfclJ  jGT ;_NUMN j GO  TO  nog 
' KPQTf t 

552  TygDE * | ) • T<M) 

_ «0Q  HcOUn*  « WCOUNT,! 


UJMCOUNT  .GT,  0)  GO  TO  553 
MCQUnT  • 40 

WK1Te«4,3P0C»  NDATE,HEU, (L IF T ( J ) , I ■* , NnODE J 
3000  FORMAT  I • i * • 2A4 , 5X# *6^4/ • g r*T  IDEST  nOM|ERI/7X|24I5|> 

553  wr J Te 16/2930 ) PAX, (TNODfeU ) , l»*»KCUT > 

If  6NP|_qT , «e i 0 ) 00  T0  554  _ _ 

NOBAf^tjODAT*! 

203o  F0BM*TiF7,2.<Z4F5.1}  ) 

554  CONTINUE  

" u«ll«* 

IF { ISPAN, GT , J) GO  TO  770 

IF  {KK»NE.NovL>GQ  TO  759 

W!flTE"(B8,780t 

780  F0BM*T»"  SPAN  1 VALUES-POSSIBLY  FCUNnATtON«) 


•0007*20 

•OOO7IM 

•0007160 

00007190 

•0007160 

•OOOtIiO 

00007100 

90007190 

00007200 

•00072*0 

00007220 

00007210 

smiiit 

000072*0 

OOOO7276 

00007210 

000072*0 

00007390 

009073*9 

00007320 

00007330 

00007360 

00007350 

000073*0 

AO6073Z0 

00007380 

000073*0 

00007490 

000074*0 

10007420 

OO0O74S0 

00007460 

00007400 

000074*0 

00007470 

00007410 

OQ007490 

00007590 

00007510 

00007520 

00007530 

00007560 

30007590 

000075*0 

00007570 

00007580 

00007585 

000075*0 

000075*1 

00007690 

•0007610 


A16 


3279T  Qi  06*11*78  lt,W 


?7i  iF|KN*ae,i»cYj;>«o  to  too 
do  7qi  T^iLO.Noae 


Ji»ntrjTTT- 

j2f IXIJ.2) 

J3f J X V 1 ' 3 ) 

JOTKUT4 } 

JPtU3*i0,44)60  TO  700 

Ttgt»H(T|jn*T|U|)^T(J0)<T(J<))44. 

GC  TO  701 

750  TSF(j)f(T)Jl)*T(j2)*T(JO))/3, 

791  IF  tNGFLG. BO* 1 ) GO  TO  752 

7»flTTf<ltT.r*r(i>-TrI.J^|Ty«eT,KCTL«TLA5STl/<D7.»OTl.> 
too  cgurutt 

Wl|TBll0.700napAM>irT>QTxPCTl.fL*CIT>P»T . BCYUKU 

760  foBHATi///*  SP*N*----— »#14,/*  SFT-HOORS-T? * FH •! * / 

1"  dt-h9u»S*-»"»fio,3|/"  sfT/oT— ,.-*,Rio,i;/ 

2"  PLACB  TBHPn"»flO,3i/"  DAYS T»El9i327 

— 5"  kdt-cocnt-.", ii/4) 

HB^TB»00I 7611 ( • #TSF,( l ) il*R|LR»NUI1E ) 

7Ql  FBBI4ATtO(tO.BlO.»>>  

UrTYE<«*,767)C!,TsF<1>,J*NELP.  NU  MF] I 

767  FBBMAT<2XM8lFlO*3) 

700  COpT ihHE  

lrUt«lWfffT 600i050,590 
090  WBiTE  16,2001)  DAYhniT<N)iN»HNUMH) 

‘ ilt.Efi.O)  - 


000076|0 

•9107610 

rh;:ii 

00007660 

00007620 

80007600 

100070*0 

00007700 

00007710 

000077JO 

00007710 

Kk;;ii 

00007760 

00007720 

00007700 

000077*0 

00007000 


V 


lFtN*i*Eli, 

m 


>00_j0  8020  

ay,in7T(N),k*1#kupn7 


WrnilVZJBTTff 

3020  LL«0 

PBAKTBMP»0 , 

DB  820  N,ji77UH|r 
lFtT^).0T,P6AKTBHP)G0  TO  026 
GO  TO  029 

■ 826  PeAKTSHP.Tr^l 
NODEpBAK'N 

P6AKDAY«DAY 

*29  Gpji/T  jiytfB 

HBlTB»6,2600|DAYiPEAKTbBP»NCDEPEAK 

c ; 

SVO  CDBTpTGF 

NBl.PfN«MB*l  ■ 

RSTURN 

C 6 • 6 • • • TO* a • *• .*•••« #S» 

c 

2001  FOBHATt/Z/OH  TJHi  ■ ,F»,2>6F  DAYS/<10 < I’* r«, * ) > ) 

2002  rofiWAT  r23H-B*7)D  TOO  L*BGE^EU7HC~, If) 

2003  F0gMATl34*OZERO  OR  NEGATIVE  AREA  ELEMENT  60.19) 

2004  FQBHaMBh  TIME  ■ .F8,2,6H  C AYS,  100X , 1H1/ < I ( 16iF7 ,1» # 35X , 1H1 1 > 
2*  OB  FWT7^r^T»«f*«PB'AK  TEMPER  ATtnr^TM  I J CAUCULATTBNY, 

tn  the  peak  temperature  cor  thet.eoui"  day  r, 


1 •*•*••••, //*  THE  PEAK 

2 r*  ?4fM."  dbqrees  at 


NODE* 


>J«) 


END 


1000701 
00007061 
00007010 
00007020 
00007060 
00007010 
00007060 
00007020 
00007? 00 
00007060 
00007920 
00007930 
00007960 
00007910 
00007960 
00007920 
00007910 
00007990 
OOOOSOBO 
00006C60 
00008020 
00008030 
00000060 
000080*0 
00008060 
00008070 
00008080 
000080*0 
00008100 
00008HO 
00008120 
00008130 


on  o 


3279T  01  0»»H-78  16.456 


i 


i 

i 

i 

i 


SQ«R0UTIN«  schsol  <kkki 
M*yNt5M»WAXWM»y#Ml»4|Q  _ . 

c 

C0HM0N  /SVMARG/  NN.MM.A(MAXK,HAXBfc);B{HAXN} 
C GO  TO  Tl OOO.tOOO  V.KKK  ~~ 

C 

C BM«6I.HA.taiI 

C 

1000  DO  200  N'liMN 

DO  260  Lji2jlmM 

C*i<H?V/i>lN»lI 
l • N4U-1 

i'  ^N*ll-jM^.ni24D 

240  J*0 

DO  250. K*U,MH 

* 3) jX < f * j| ”P • A ( N»  K)  ' 

269  A»N.C)«C 

980  CoMTiagE : 

GO  TQ  500 


RfePUCB  VECTOR 


3009  Db  200  N»1*MN 
08  205  L.l.HH 


I *N*U,1 

IFJNN-I)  290*285 j 285 
295  iUi*0  II 1 - A<H-»ia«J.<-ia- 

c *’•  S.TUMStWtt, 

c 

N«NN 

300  N * N^i 

__  __  LEL1NJ  35  a,  SO0.35J5 

TWO  DO  4gO  k«2,hh 

L » N4K-1 


»tN>  « 

400  CONTINUE 
(30  TO  JO 0 

500  BBTURR 


Bill 


END 

SUUROUTJue  TRLOT 

_BftUBJL_  

end 

BfiECUTE  0UM6 
blHtTj  40;33Kt;25000 
E tUE  07,Dlg;iu 


EJUE 


98,829,61 


0000*140 

99108108 

10808180 

00808160 

00008170 

00008180 

80000160 

1MIHI 
00008220 
_0 00 00? 30 

80808210 

89008210 

—88808269 


00008270 

0000828b 

00006?*8 

80008300 

00008310 

"JS88SH! 

00008340 


ooooejig 

00008340 

00008370 

80006380 

000083*0 

000064QC 

00008410 

00008420 

00008420 

00008440 

00008460 

00008460 

00008470 

00008480 

000084*0 

00008599 

000085 1C 
00008520 
00008538 
00908540 
0O0CB5I0 

00908560 

00008570 

00008580 

000085*0 

00008690 

00008610 

00008620 

110008630 

00008640 


I 


I 


A18 


! 


I 


APPENDIX  B:  TWO-DIMENSIONAL  THERMAL  STRESS/STRAIN  CALCULATION  PROGRAM 

(WES  VERSION) 

I 


i 


B1 


n n no  noon rii 


1UST  01  04*11-74  - 15.4lZ- 


ident 

MSG  3 
OPTION 

forty 

LIMITS 
FILE 


R0CC46, CAMPBELL 

062678/1700 

fortran 

XKEF.MAP.PEGK 
,33K, ,6000 


C». X3S, 1L  < NtWrBSTBESDW 

STREsSV-WITH  VARIABLE  tjme  and  strain  cap  input 
••••copy  AS  OF  27  MARCH  704NCL  STRAIN  CAP,*  SPECIFIC  ANALYSIS  TI5IS 

ARBITRARY  twO-DlMENSlBN*L  STRESS  STRUCTURE  INCLUDING  INCREMENTAL 
CONSTRUCTION,  MCHfcNHY  CREEP.  RESIDUAL  STRESSES,  THERMAL  STRESSBS , OGOgOlSO 


OOOoOOlo 

00449015 

OOOoOOlo 

0000°ORo 

OOOoOOlo 

OOloOOlo 

OOOoOOJo 

OOOoOOlo 

OOOoOOlo 


VARYING  PRESSURE  BOUNDARY  CONDITIONS,  AND  BlMODULAR  HATERJAL 
properties. 

HXN»HAX  1MUM  _ MQJ__QL_MttCES 

MXE*MAXIMUMUM  NO.  of  elements 
MN2*hAX  NODES  « 2 


PARAMETER  MXNs576»Mxt«475,MN2*MXN*2 
CHARACTER  NX7*6(7),N77*6(7) 


OCOQOlfQ 

ocoooiio 

OOOoOlIp 

OCOoOlli 

00000182 

OCO0OIB3 

OCOoOlfA 

ocoooiBe 

OCOoOlfo 

OCOflOllo 


1, NgT,NCClJNT,TI,DT»DTT,TliTL,Xi:,YC,ST<3»10),TIMLAt27).NUME(27>.Tl(.,OCOo016o 

2TTT(15)  .NUMOL, TIME.NNAL.NL Ay, I SC, EX, EY, DASH, ANAL (200>.NANALT,T LAS  TOCO  0 017o 
3.RRR  (5 1 , ZZZ15J  j NFLftGsILl.  — - - OCfioOlBo 

COMMON  /MaTaRG/  E<3P,6.8)»RO(8>,EB(I),HEDI15),CICI20,4,»),  OOOoOllo 

1CC(4,8),NCREEP(8),SCL30,2,9),NSC  OOOo02lo 

CQMMQN  /ELEARG/  IX  iQxE . I ) xMTAGiMxB  )IS  IQIMXEt?  >-,TQLD  tMXE ) i 5004021c 

1DE1KMXE) ,DE12(MxE),DE21<HXE) ,Dfc22(MXE>.DSIGT6>,CC0<4>,CCC(4?  000fl°2fo 

2, EES<MXE,2).EPS(6).SCAP(MXE)rNE7(7>.NX7,NY7,YEMPE(HxG)  OCOo028o 

COMMON  _/PR_SARG/..lBCAiT!01.i.JBCI  lflP)_rPR  CtOS)  _ _ OCOo02|0 

COMMON  /OHDARG/  R(MXN),7(MXN),LR(MXN)*U2(MXN),C0D6(MXN)»T(MXN)  00000200 
COMMON  /BANARG/  MBaND,NUm8lk ,B(120)Ia(120#60>  OCOo02Bo 

_ COMMON  /l  S4ARG/  1 ,J.K . S < 1 fl  ■ 1 (n  C J 3. 3 ) . 013,_31iB  1 5 . 3 ).P  (.10  )^tt  1 4.U  044flQ2lc 

1 F ( 3 . 3 1 0000  02|o 

CHARACTER  REaDM0DE#6«DASH*6 , 0VER«6  000o02*o 

SFTIM*8,  00000305 

dash»6H  comp  oceooaio 

QvgR«.6.H_>lKXi , 99440314- 

NF«20  000o03|o 

WRITE(7,4100>  0000°3»0 

RFAPHoDEiAH  0444031c 

c«*«**read  and  print  of  control  information  and  mayereial  eRoRV  oooocsji 

000fl03|2 

SO  READ<Nf.l000.END»5g0P>HEB 0040 0.314 

OCOoOSJo 
00008300 
005ij03f  q_ 
OOOoO40q 
000o04|n 
000q04|o 
OOOo048o 
OCOq04|o 


C.»READ  MODE  of  INPUT-  "TaPE"  or  "CARD"  FOR  NODES, ELEMENTS.  4TEMPS. 
RFAD(NF,998)READM0DE»NEXTRA 

___  READ  ( NF  , 999  >NuMNP,  NuMEL,  NUMHaT  ,NLaT  ,ISC,NaNALT,TI,Ti 

WR i TE ( 6 , 1009) 

WRITE  {672000)  HED, NUMNP, NUMEL, NUMMAT.NLAY, JSC, TI.NaNALT 
56  DO  59  M.l.HUMMAT 

REaD<NF,1001)  MTYPE,NUMTC.NSC,NCReEPIMTYPE>,RO(MTYPE) 


B2 


k 


-341**  -at-  01-11-2*  15.617 


WR j Tfi  (672011)  MTYPE*NUMlCi.NCRE6P<8T»PE>jR0(nTYFE) 
WRITE (6, 2051) DSC 


READ (NF. 1005)  ( (E ( 1. J,HT7PE>.J*1,6>. W.NUMTC) 
WRITE  (6.2010)  (<E*I»JiMTyPE).J*1#6)II*1> NUMTC ) 

00-5*  liWUHICj-SQ 

00  5*  d»l,6 
5*  E( J«HTYPE)*6<NUMTc»J#HTyPE) 


RE  AD (NF,  10  08) ( <SC< I. J.MTYPE)7j«l»2).l»l.NSC) 
WR!TE(6,2052)( (SC( I, J.MTXPET;j*li2>.l»l.N8C) 

D.Q.  75  l.«NSC*30  ._.  . . 

DO  75  J«l . 2 

75  SCIIi'J*HTYPE>»SC<NSC«JiHTyPE) 


IF  (NCREEPIMTYPE))  58.59i54 
54  NCR«NC8EEP(MTyPE> 

RE_.lDINfjJ.QDai-  ,aJTlUt(Cl£lUJJuPMPE)ia?1.4)»I»liJllCRl 

WRITE  (6,2013)  HTYPE.  (TTTU),(CIcUIJ*HTyPEJ,J»1,4>.I»i*NCR> 
REaDINF,i005>  <CC<l,BTYP6uipl>4) 

WRITE  (67201*1  t CC(  I •HTYPE) . Jtlj4) 

59  CONTINUE 

c«**«»read  and  print  nodal  loaP/D I splacehent  BOUNDARY  COND?*#* 

WRj.TLl_6j2.QMl .... 

L»i 

3006  READ  (NCI  3000  IN,  COOjjl  N)  »U*CN)*UZCN) 

lriN.EQ.DGO  TO  3006 

3005  L»L*1 

IFfN-Ll3002.3001.30o4 

3004  C0nE(L)«C0DE(L-l) 

UR(L>»°. 

UZfllsO, 

on  TO  3085 : 


3002  WR!”E(6,2017)N 

call  exit 

IS 01  I F (L , EQ . NUHNP ) GO  T°  3006 

GO  TO  3006 

C«**«*READ  NODAL  POINT  DAT*  FROM  TAPE/DlSK  AND  PRINT* 

. aoo*  L»i_ 


1F(REaDH0DE,E0,6hTAPE  )fiO  to  905 
READ (NF,  4000  HM, R (N ).Z (N)  ,N*1,NLMnP) 

fiQ._-IQ.-iM , . 

905  READ(4)(m,R(N5.Z(N>.N«1,NUMNP) 

906  CONTINUE 

WRITEt6.2002)(K,  COPEtKl . BlK ) »2  (KI-UR4K)  , UZ(K ) »Kll«  MuHNP) 

WR!TE(3)(K.CODE(K)«R(K),Z(K).UR|K)»UZ(K)>K«1tNUHNPI 
WRITE  (6.2001) 

C*»«_«#READ  L ELEMENT  HA  IA  ffiflM  TAEL  QB_  CARDS  AND  PRINT  

NmO 

IFl REaDMODE ,EO. 6Wr *Pfe  >?0  TO  90* 

. READ (NF. 4001) (M, (1A(N. I)jl»it5). t SlQlN, 1 > . I il . 3)« N«i»NUH6L ) 

GO  TO  909 

908  READ(4)(N,«IX(M.I ) » I *1 .5 ) » (S IG (P i 1 1 . 1 «1. 3 ) , NAIiNUMEl i 


B3 


00080450 

ncflfl04*p 

00000470 

OC00046o 
00000490 
OOQoOSf 0 
OOOoOSto 
00*6  03  to 
00  Oo  05lo 
000005*0 

00000550 
000005*0 
00000570 
DQ**05M 
00000590 
0000  06*0 
oooooof 0 
0C00°8l0 
00fl0°8|0 
00*0*6*0- 
OCOo065o 
000006*2 
OOfloOOBo 

0000°67o 

000o06*o 

OCQQ0690 

ocoo076o 

00000710 
OCQ*  0 7*o 
000007*0 
000007*0 

_D0  Do  0-7*0- 
000o07«o 
00000770 
OC DO 07fio 
000007*0 

0C0007»i 
OOOnO 8*o 
0 0 00  0 8*0 
OC  00  082 o 
0.00058*0 
000008*0 
000008*0 
OOOq  08*ft 
006oOB7o 
00008860 
0.00086*1 
OOOoOSfo 
0000091 o 
Q0Dn09ifl 


00000920 

000009*0 


lllST  01  OO-li-TA 


909  CONTINUE 

HR l Tf ( 6 , 2 0 03 > ( N , ( JIU * 1 flumWILL I It  I •! • 5 Lf 

HrTtE < 3 ) ( N . « Tx I N , IT, rflTTuTsTGlN;  I )II»  1 • 3 > # N»  1 , NUHEt ) 
j»6 

00  3<0  N*1 , NUmEL 

00000900 

0000°9T0 

ooj(jove0 

DO  340  1-1.4 

Do  325  L-1.4 

KK-IABS< IX«N, I »-IK«N-L>  > 

00000900 
OCOOlOlo 
OOJolOfo 

IF  (KK-J)  325,3257^20 

320  JsKK 

325  CQNT I NUF 

OOO0IOE0 
OOO0IOS0 
. OQOolOPo 

340  CONTINUE 

HBAND«2«j*2 

nunol*o 

00001050 

OOOOlOfo 

OOOOlofo 

nanal«o 

TIME»Tl 

tlast.ti 

OOOOlOfo 

OCflolOTo 

OCB0UI0 

ocoollfo 

OCflollf 2 
OOflfllllo 

3026  REaD(NF,302C)M.TolD(N) 

IFIM.EO.DGO  TO  3026 

3Q25  L«U*1  

OCOollOfl 

ocoouio 

ocooiifo 

ocoollfo 

OCOflllTo 
....  OCfiollOo 

IF(M-L)3022. 3021. 3025  ' 

3024  T0LDIL)»T0LD(L-D 

GO  TO  3025  

3022  URt TE(6.3027)M 

call  exit 

1021  IFcL.EO.numEUGQ  TP. 3o3i  _ ...  

GO  TO  3026 

3031  L»0 

HR  I TE  < * i302« ) ( M,_iOI,DiH  > ifl-J , NUBEL ) 

OOOOllfO 

00001210 

OC0012IO 

OC0012»0 

ocoolzfo 

00#ol2fo 

DO  450  Nrl , NUMEL 

DEll(N)«n. 

DEi2(N)«0.  

DE21 ( N) »0 . 

DE22IN) »0 . 

SlG(N,4}rO,  ....... 

OOS01250 

0C001260 

0Cfl0l2?o 

00001280 

0C001290 

OC0013«0 

SlO<  N , 5 ) *0 . 

MTAG(N)»1 

0C0013I0 

OC001320 

I F <SlG(N,l>)  445.445.44P 

440  IF  (SlG(N.l)*Sir.(N»2>)  441,441,442 

441  HTlG ( N ) *2 

CO  TO  45 q 

442  HTaG(N)*3 

445  If  ( S iG (N, 2 ) ) 450.450,448 

448  MTAG1 N ) *4  . _ _ . 

450  CONTINUE 

00  4*0  N.l.NUMNP 

Fri2»N-i)«o. _ 


OC001S80 

ocoolsio 

OCOolJSo 

OCOOl3|o 

OCOOl3T0 

ocooisfto 

OCflOl3»0 

OCOOl4#o 

00001410 

OC001480 


460  FFf2«N)*0. 

C*****REaD  ANALYSIS  TIMFS(0AYS>  FOR  RtN 


OCOOl4|o 

OC0O14I1 


B4 


JU8T  Ol  aa=U -74  15.6 1? 


READ (NF, 3017) < ANAL < N > » N«l,8ANALT > 

HmTil6i3042)  ( ANiL<N*  .N»1I8aNAL1J 

3042  F ORMAT f ///"  ANALYSIS  TIMES  FOR  RUn(DAYS) ■// , llX, 18F7?1 1 ) 

NAL»1 

NCP*Q  — . 

IF(NExTRA,E0.1)NCP*1 
DO  8?0  KK-l.NUMMAT 

IFtNCREEP(KK) .GT.fl)MlP»l 

#70  CONTINUE 
890  DTT»1.0 

NfL  AG»0  _ _ _ 

NANAl»l 

C»«»»#READ  LIFT  CONTROL  

DO-  800  -LLLalxNL-AY 

550  READ ( Nf > 1006 ) L AY, NuMN. NUmE ( LAY  1 ,NUMPC*NP. NDT . T IHLAcLAY ) *T IMNL 
WRITE  {6,2008  ) LAY. NUMN, NuMEC LAY > JNUMPC, NP, NOT , T IMlACLAY ) , T I0NL 

JFLAGiO _ - _ 

NNALai 

nunnl»nume(lay) 

IF(NFLAGjQ.2)aQ  T9  ?51 

897  IF(NNAL,FQ.NDT.AND.NCP.E8.1>NruAG«l 
C«»«**REaD  NODAL  TEMPERATURES  FrQm  TAPE  OR  C6RDS»*»»»ai»a»*4*t» 

’5i  IF (RFAPM0DE.EQ.6HTApB  )Go  TO  901  _ 

READ(NF.3009)TIMETAp8 
READ (NF. 3007) ( T ( kK) >RK«l«NUHN) 

GO.  TO.  .?  02. 

901  READ (4*3040) TIME TAPE 
READ(4.304l)(T(KK)«KKsl,NUMN» 

902  CONTINUE  

IF*  NAN AL I GT . 1 ) GO  T0T40 

. IF<TIME-ANAL<NAL))944,948»?43 

*44  DT»TImB-0. 0 __  _ 

IFCTIHBTaPE-TiME)94i«952*941 
«40  IF<NFLAG.E0.2)G0  TO  947 

IE (TIHETAPEiTI ME >9^*952*953  

947  IF (TJMfeTAPE-TIMETJ 951, 952.953 


952  DO  871  KK*1.NUMNL 

I*ixlK£»JJ 

J«IX<KR»2> 


K*J  X ( KK . 3 ) 

L*IX1KK«4> 

irc'K.ED.UOO  TO  872 
T£MPE(KK)»<T< I)*T( J)*TIK)*T<l))/4 

GO.  

872  TEMPE(KK)3(T< I)*T(J)*T(K» )/3* 

871  CONTINUE 

IF(NFLAG.NE.2)G0  TQ  095 
DO  874  KK*1,NUMNL 
874  TSAV6(KK).TEMPEtKK» 

895  IF(NCP?NE.l>GO  TO  95? 

975  DDTl»ANAL<NAL  + l>-ANAL(NAl.) 
IF('nFlAG.EO.O)GO  to  596 


OOOoHlo 

DC4fll4fo 


OOOoKBo 

00001470 


008ol4!i 

00801472 

00801413 

DOOo14Ia 

0C001475 

00801476 

00801477 

00801478 


00001479 

Ofl|0l48o 


00001490 
0000158 0 


0C00l5ll 

00001515 

008015*0 

.00001515 

OCOol52o 


00801585 

00001550 

OCOol58o 

008ol5?0 

00901580 

000015*0 

0080l&8o 

00001810 

OOOq168o 


ocg0i6|o 

008015*0 

ocooioJo 


00901855 


OOOq168q 

ooeoioBi 

00001682 

00001083 

OC00lO6< 

OCO0IOB5 

00001086 

OCO0IOI7 

ooooioie 

.OOOOlOIfl 


00001672 

00001074 

00001675 

00001676 

0000lO?7 

ooooiora 


ooooioio 

OCflolOll 


mar  01  08211-71 


t 


GO  TO  {960  <985, 964*901*  903) /N0LA0 

961  DO  014  KK»1.NUHNL  . , 

00001802 

00881804 

804  TEMPE(KK)«TSAVE(KK) 

896  IF (DDT1 , GE , 1 ,0 • AND* DOT 1 , CT,2)0TT«0*5 
IF(DDTl,QE.2.)DTT»l,0 

000*18(6 

000016*0 

GO  TO  *62 

905  NEw»N|jME  (LAY-1) 

DO  880  KK«1,NEW  „ _ . . 

BOO  TEMPfTRKV*(TSAVE(Kt<)»TOL0(RKn*0rTXN?24.tTOLD<KRr^ 

. DO  881  KK*NEM*l>NUHNh 

881  TEMPE(KK)iTOLD(KK) 

00101808 

SliSiSU 

GO  TO  *60 

964  DO  883  Kg«l*N()HNL 

®83  TEMPI (KK  >»(TSAVE(KK)»TOLIiRRT)» f5«TOtBtKK>  . ...  . 

00101703 

962  CONTINUE 

lF<Nf.  LAG  .BE.  3i50-JtQ_Ufi 

If »-■  — - 

955  WR 


iriNFLAGlEO.O, AND. JfUO. 6091100  To  3(0 

WNITE(6, 2009) NANaLi UlTi TlHCTlPEf (KK, } < KK) f KKiljNUNN) 


C»*««»READ  P8ESSURE  b.  c.  From  cord  or  t«re*«< 
290  HRfTE  {6,2005) 

DO  300  UlJillMPC _ 

READ  INF. 1004)  !BC(t )«JBC<L> , PR(l) 

300  WRITE  {6/2007 J 1 8C<L > , J80a>  »RR<U 

3 10  CQKTlMUE 

JFLAQ.O 

WRiTE(0,«009)NFUAGrDTT 

4009  FQBM*I{/"  NEUG«»al4s"  _ 

!F{NFLAG.GT,1)G0  T®  55 
WRITE ( 6, 2025 ) 

2025  FQHHAT{/;»  •pflODAL  PtSPLRCBMBNTSVl ) 


>*<»• 


55  CONTINUE 

IF  (NPrl)  435,435,500 

435  00,350  J^liNUMNI - 

350  MTAG(Nl«l 
500  CONTINUE 

NCOUNT80  : 

DO  570  NNN«1»NP 
425  NCOUNTSN COUNT »1 
C«*f*F0RM  STIFFNESS  T1 
CAUL  STIFF 

C»««**CalCUL»TE  NODAL  DlSPLACEMEBTSA«9*«6«i*««»***S«»«W#*#?**« 

CALL  BAhSOL __ 

IF  (NCOUNT-NP)  525»510»510 
510  DO  520  N*1 » NUMN 

C»*»«*IN!T!ALIzE  DISPLACEM6NTg4»*.s**«#0*M*t*«*«*8«««»*»«****  . 

FF  ( 2»N-1  > *FF  < 2»N-1~>  *8  ( 2*N-t ) 

520  FF(2»N)*FF(2.N)*B<2»N) 


IFtNEuA 
WRITE  { 


GiGT.DGQ^O  $25 

6,2006)  (N,EFW*N*lI.rr(2*N).H»l#NUHNJ 
WRITS! 3) (N.FF(2*N-l)«FFt2«N) /N=1,nUMN) 


WV0O*'»7 

OC00l7|S 

00001719 

001017(0 
00001702 
_ O0#al7to 

00001705 

00001700 

00001700 

00l0l7fo 

00001700 

.iifloai7fs 

00001701 

00001705 

00001706 

00001707 

00801700 

_oejfti70D 

000017*5 

OOOQlOlo 

OOflOlOfc 

00001800 

00001800 

_OGOalBfO- 

000018*0 

00001800 

00001801 

00001800 

00001801 

00001800 

000019*0 

0000190b 

00001901 

00001900 

00001900 

_00Jo19J5 

OCOol9Ro 

000019*0 


L 


V 


suai  01  o*3u-7a.-  15.61Z  — . . 

994  Cai  l CiEEP 

OQOOltfi 

OOf glVlg 

370  continue 

NU»lOL«4UtlE(tAy)  . 

IF((NaL»1).QT .NAN ALT >G0  10  9000 

jr{'Nrt*G;QT,oyNFL*«»nru»*t 

ir (NFLAfi .Efi.41 JFLifial 

OOOol’IO 

ooigiah 

mum 

001019)4 

juutuU 

I F ( NFlAG . EO • 0 ) GO  to  969 

GO  TO<967,»67.*67.96?.97B>iNFLA0 

•6*  NfLlGaO  - — 

Dt«TLA8T*1.-TIME 

tim«time*dt/2. 

time«tlast*i. 

00001916 

000019*7 

00001900 

00001902 

00001904 

00081906 

TLAST«IIME 

DO  971  KK«1.NUMNL 

971  TFMPB<KK)«TSAVE(KK>  

GO  TO  900 

967  Dt»SPTlP/24, 

Ttm«T!ME*DT/2. 

00001908 

OOOOZOIO 

00002009 

OC0Q2O0O 

00002009 

00002000 

T!MET«TimE*ANAL<NAL*1)-ANAUNAL> 

TIHE«TlME*DT 

1F(NF1.AG.E'1.2)GQ  TP  098  . ..  

00002002 

00002009 

00002007 

GO  TO  979 

969  !>T« ANAL <NAL«l?-ANAl<  NAD 

tIM«TiHE»QT/2. _ ....  _ . . 

OOOO2O0Q 
00002009 
00002019 

TimEaTTME-DT 

TLASTaTlME 

898  NAI,  «NaL*1 

00002090 

00002093 

00002009 

nnal*nNal»i 

nanal«nanal«i 

400  IF(NNAL.OT.NDT>GO  To  600 

00002080 

000o20>i 

08002002 

IF(NNAL.EQ.1>G0  to  899 

GO  TO  897 

600  CgigTlNUE 

00002016 

00002088 

OCOo20To 

00002000 

00002090 

00002100 

GO  TOTO 

94i  WRtTE(6i4002lTIME*Tl tlETAPE 

4002  FORNATt"  ABORT— FIRST  ANALTSIS  TIm6  AND  lAPETIME  DO  NOT  "1 

" 1" AGREE.'  »/»  Ting—-.  > -aF7;2/-  TICETaPE  * SF7,2> 

943  WRTTE(6,4003)TIME.ANAL(NAt) 

4003  FORHATC"  ABORT— TIMES  OF  FIRST  ANALYSIS  ?IVEN  BY  Tl  AND  "i 

0C002110 

00002120 

OOQo218o 

1"ANAL(1>  DO  NOT  AGRE6:"/?  Ti*r-— . ■ «,F7.2,»  DAYS"/ 

2"ANAL(1>-  * «;F7.2."  D AltS'J 

Gn  to  Sono 

00002190 

00002190 

00002100 

953  WRiTE(6,4004)TlME.TlMETAPE 

4004  FORMATC"  ABORT— next  ANALYSIS  TIME  CANNOT  BE  FOUND  ON  TAPE", 
i/»  TIME-.—  « " ,F7 . 2»"  DAYS"/"  JImEIAPE.  T "iF7,2/"  _DAYS"1  - 

OC9o217o 

00002180 

80002190 

GO  TO  9000 

963  WRTTB(6,4006) 

4006  FOBMATf " NFLAG  CanNoT  ■_ 5 _8T.  SlATEflEMT  9Z3*"L  _ 

00002200 
00002201 
. 00002202 

GO  TO  5000 

970  WRtTE(6,4007) 

00002203 

OOOQ2204 

JU8T  01  DI:llrM 15. 


4007  FORMATt"  NFLAG  CANNOT  » ( AT  STATEMENT  9#9-») 

' 998  ^RHAt!?“a17T51 
999  FORHATC615.2riO.O) 

1000  Fo«HATtl5A4)  

1001  Format {4 i 5»Pio .6 > 

1002  FORMAT  ( I5.F5 . 0 » 4Fl0 • 0 ) 

1003  FORMAT  trig. Q.4E1Q A 

iooi  format t<Ji5t fi o .0 > 

1005  FORMATi6FlQ,0) 

1006  format C6 i 5 . 2Fio^o) 

1008  F0RMAT12F 1 0 . 0 ) 

1009  format iiHi > 

2000  FORMAT  (15A4Z 

1 40HO  NUMBER  of  NODAL  POlNTS-r-»-.r-«— 

2 40HO  NUMBER  of  PLANTS— 

- _3  .4J0H0  ..NLMBlR_QF_n.lFjUJiEM!_ilAl&KlALS-k--— 

4 40H0  NUMBER  OF  LAYERS  IN  THE  STRUCTURE  — 


5 4 OHO  STRAIN  CAPACITY  INRUT 7 1-VE§.  8-NO  — 15/ 

6 40HO  TIME  OF  FIRST  8NALrsIS—-a^— TIBIA/ 

7 40H0  TOTAL  NO.  OF  ANALYSES  REQUESTED— I—  14) 

2001  FORMAT  (92H1ELEMENT  No.  I J K L MATE 

1G1-RESIEUAL  SIG2-«ESIPU6L  ANGLE  ) 

2002  FORMAT  ( U2,Fl2.2.2Fl2.3i2E2417) 

2003  FORMAT  (1113*416. lIl2<2F17.3«F9|3) 

2004_F_0RM.AT_.  (97H1N0DAL  POINT  TYPE  ..XrftSOlMAlE-  Y.QRdJNATE 

1AD  OR  DISPLACEMENT  T LOAD  OR  DISPLACEMENT  ) 

2005  FORMAT  ( 79H0PRESSURE  BOUNDARY  CONDITIONS/  24H  I 

1URE  > 

2006  format {5 1 26H  n.p.  displ-Ux  display  >/?< i5*2Fio.5,ix) ) 

2007  FORMAT  (2I6.F12.3) 


MATERIAL 


2008  Format  (SQHi  NUMBER  Of  LAYERS  IK  THE  ANALYSIS— ‘—-..C.-r— *15/ 

1 50H0  number  of  NO&AL  POINTS  IN  ThE  ANALYSIS— —6— -=«  15/ 

2 50HO  Number  OF  ELEMENTS  IN  THE  ANALXSISr A--_-«i5/ 

3 50HO  NUMBEB..Q.F  .PfiESlURJE^CARfiS  ££R.  TH£.ANALYS!S-i--— *15/ 

4 50H0  NUMBER  OF  APPROXIMATIONS  FOR  STRESS  CALCULaTIoN.,15/ 


6 5oho  number  of  analyses  at  this  sTase  of  constr'— -»«i5/ 


00002215 

-Minus 

000022(0 

000022(0 

OOOq228o 

000022(0 

— mm 
000022(0 
00002Z?0 

00(023(0 

00(02310 

008025(8 

000o23(o 

OCflo23(o 

OCQC2350 

OGOo23(o 

OCBo237o 

0 0882  388. 

00002390 

.S!OCOo24Sg 
OC0024IO 
OCflo2«Co 
00002«(0 
X .LOOC8o2Alo 

OOOo24(o 
PRBSSOCOo24>o 
OC0C24FO 
OCOo24(o 
0000249c 
15/  .ocoo25(a 
0C0025J0 
00002520 
_ 0000258c 
OOOo25Ro 
000025$0 


7 5QMCL  .1IME  QFLLAYl.Nfi._JJl£  JTflp._UUar-»— z-— — — 6-j1jrE10*3/__  . 0C0o25Bo 

8 50H0  TIME  OF  laying  THE  NEXT  LIFT— r— * — --.«P10,3>  OCOfl25To 

2009  FORMAT  (42H1  NODAL  TEMPERATURES  FOR  ANALYSIS  NUMBER  15.  0C002580 

1 " STRUCTURE_UP  to  L TFT  ?,I5«"_  TJME  FRQM  tape  «-,Far2*  _ 000025(0 

2 » DAYS"//"  NP . TE0P",?(13H  NP-  TEMP)//  OCOo26»o 

3 10 ( I 5*F8 , 2 ) ) 000026*0 

2010 . format  u shq  temp./timE  iqx  sheicj 9x_«lhnb Ux . aheuj o too 2 6 g a 

1 ldX  5HG/H2  lOX  5HALPHA/  000o268o 

2 (Fl5.3,4Fl5.5.El5.5>)  000026*0 

2011  FORMAT  (i7H0M>T£bla^  NUMBER*  13«  SJBHl  NUMBER  CF  TEMp7/TIME  CARDS  «OCBo265o 

1 1 3 . 24H*  NUMBER  OF  CREEP  CaRDS*13.  15H,  MASS  DENSITY.  El2;4)  OCOo26(0 

2012  FORMAT  (26HONODAL  POINT  CaBD  ERROR  Nk  15)  OCOc26To 

20J.3  F.Q8BA.T.a7J18_MAT.E8lAk_NUMB£R..I5Z/ 0C0o26Bo 

1111H  COEFFICIENT  FUNCTIONS  A < T ) IN  MCHENRYS  ECUATIOn  STRAINJI)  « 000026fo 
2STRA JNiO>*Al(T> (l-EXP(-Ml»TJ)«A2<T) (1-EXP (-M2«T>)//  OCOo27iO 


;_4ME<t>  _ 


w 


Sim  01  08:11-70  15.6*1 

310*.  10HTEMP./TIME  llx.  24UA1.A2  FOR  COMPRESS. CREEP  12X.  23K0I.A4  O00027f0 

4F08  TENSILE  CREEP//  3flXtJJ2ttAlLljx.2HA2. llx. 2HAS.I3k:2HA4// OMoEJlt. 

5(iaX.FlO.3.1OX.ElO.3*5X,Kl0.3,lOX;ElO.3.5X.ElO.ST)  OCOg27Bo 

2014  FORMAT  (30H0  INDEXES  IN  MCMENRYS  EQUATION//  OOOO^iO 

13QW  For  COMPRESS l Vt  CREEP  Ml  «E10 . 3«4X » 4Hh2  a Eio?3/  00002750 

23 OH  FOR  TENSILE  CREEP  M3  ■ 610 . 3, «x. 4HM4  * BlO.S)  000027lo 


2017  FORMATt"  ERROR  IN  Pl6P/0«  FORCE  SpBC.  AT  NODE-  *JIS> 

2051  Format 0*  no.  of  stBaAn  capacity  PtS.  » «.M> 

2052  FORMATt"  TIME-DAYS  STRAIN  CAP WIN . / IN , "//2 (2X» FlO i 0 ) I 
3000  FoRMATf I5.F5.0.2F10.B) 

3007  FQRMATC12F5.0) . 

30i?  format ti2F5 , o 5 
3009  FORMATtrlO.O* 

„3Q*1— FORMAT t21 5. 3Flflx21 

3020  F0RMATt2x, I6.F10.0I 

3027  FORMAT!//-  ERROR  ON  STRESS  FREE  Te«P;  INPUT  AT  ElE , iMa  *•  * 1 5 ) 


OO0e27»O 


000027*0 

OOOO260O 

000026*0 

00002#*0 

oooozoto 

OOOq2«Bo 


3028  FORMAT!//" 
110(15. FO 1 2 > > 
3040  FORMATS, F8.2) 


element  stress  free  TEMPERATURES-;// 


3029  F ORMAT 1 1 X , a8 ) 

3030  FoflMATiA8.A8.Al) 

4000  FORMAT;  IiQ.FlQ^Q.FiB^I - 

4001  FORMAT;J5,5I5,3F10.0> 

4100  F0RMATJ20X. "PERCENT  MAXIMUM  TENSIlB  STRAIN  WITH  RESPECT  TO  "a 

1? STRAIN  CAPACITY". /2PX.73HH»>//25X.*4  - STRAIN  OF  STRAIN". 

2"  CAPACITY  > OR  s 75X" » /24X . "«•  - STRAIN  OF  STRAIN  CAPACITY" 

3"  > OR  * 90%"// ) 

5000  STOP 
End 

subroutine  quad 


PARAMETER  MXN.576.MXg. 47.5 . 

character  nx7*a<7>.n*7«6<7) 

COMMON  NilMNP.NUMEL»NUMPCiN,;VOL.TEMR,MTYPE,LAY.NUMNtNANAL«NP,OVER 
l.NDT.NCOUNT.TI .DT.DtT,T1»TL, XC,YC»ST(3.10), TTKLAT27,,NUME (27J.TIF.  0 0003056 

2TTT <151 » NUHOL;tlME*'NRAL»'SLOY,ISC*E3,EY.DASH,  aNAU20o>. NANALT *TLAS TOO Oo3OB0 


00002000 

00002800 

00002870 

OOBo20Bo 

00802800 

_B0tti291o_ 

00002910 

000fl298o 

_OC0029BO 

000029*0 

00002950 

_MAft29flO. 

00002970 

00fl0298o 

0.06fl29fa 

OCOO^OBO 

000o30to 

OC0Q3O12 

000030*5 

000o30to 


3.RRR(5)i2ZZ(5).NFL6Q»TL1 

COMMON  /MATARG/  E(30.6.8>aRO(8),eE(»>.HEP(15).CIC<20,4,B>. 
ICC  < 4, 8 i iNCRE  F>  ( 8 j V SCL30.2.6) »NSC 

COMMON  /ELEARG/  IX(MXfc,5),MTAG(Mx6):SlQ<MXE»5).T0LD(MXE). 
lQEll(MXEA,DEl2(MXE.>.PE2l<MSLl.Di22!MXE)JPSIGTJ)iCCQ(|)xCCC(4I 
i).SCAP(Mxt;)»NE7(7),NX7.NY7,TEMPE(MxE) 


00003000 
00013000 
OCQo307o 
000o306o 

0G00309q 

2;EBS<MXE72).EPS(6).SCAP(MXE)»NE^(7>.NX7.NY7,TEMPE(MXE)  OCOO310O 

COMMON  /PRSARG/  IBCUOO)»JBC(100>.PR<10Q)  OOOo31Io 

COMMON_/CHDARG/__R(M^(0.j-ilMXN)iIiRlMXNJ_tk?!MXN ) iCODfilMXNl  *7_IMJ(N  ) 00flo318o 

COMMON  /SANARG/  HBaND.NUMBLK,B(120).A<120.60)  0C003l?0 

COMMON  /LS4ARG/  I . J*K . S! 10 , 10  ) . C < 3 , 3 ) . D 1 3. 3 1 . H ( 3. 3 ).P CIO ) ,LM 1 4 > , OC0031»0 
1 F < 3 > 3 ) 00003100 

1«IX(N*1>  OCOo3l8o 

Ja i X ( N*2 ) 00003170 

K«iXtN*3)  .... . OC003180 

L r I X ! N./4  ) OC0o3l»O 

I X(  N.5I  * I ABS<  IX<  N » 5)  ) 00003*00 
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SHOT  01  0fl"ll*7.8 15. til 

HTYPE* IX ( N , 5 ) 

5***#o.FOtB  STRESS-STRAIN  BE 
C###** TEMPERATURE  DEPENDEnC 
DO  15  N«1.4 

MM»1X(N,m>  

RRR (M  JXR (HH) 


MSHIP  InCLODINQ  TIhE  QB**#T»*« 
Elastic  constants*#*#****.*#** 


RRR (M  j*R (MM> 

15  ZZ7»M>SZ<HM>  ( 

TrM t TlfH5o! 40l50LD  N >/2i 
40  DO  103  H*2 » 30 

IF  (E(M.l.MTYPE)-TEM) lOS.lOA.J^ __  _ 

103  CONTINUE 

104  RAT  10*0 ■ 0 

DEN  * E 1 M t ij  MI  vpE_i  t£AM*1j-1iJ1I  YPEl 

IF  (DEN)  70,71.70 

70  RATIO*. (TEH-E(M-1,1,HTVPE) I/DEN 

GO  JO  71 — - - 

50  DO  55  M.2,30 

IF  (E(M,i,MTYPE)-TL)  55,00*60 
55. CONTINUE 

60  R AT  1 0*0  , 

DEN«E(M,i,MTYPE>-F(M*l»l*HTYPE) 

IF.(DEN)  64,71,64 
64  R*TI0»(TL-E(M-1,1,MtXPE)I/DEN 

71  DO  105  KK*1 » 5 

105  EE(KK15E(M-l,KK*l,MT!‘PD*RAU0#c£iB.KK*liMTYJ>El-?EtN-l*XK4l*P.TTPEJ 
DO  300  KK*1 , 2 

300  EES(N,KK)*EE(KK) 

1 F l ISC- EQ. 01(50  To  76 
DO  455  HM*2, 30 

IF  ('SC(HH.I.MTyPE)-TLI  >455,460, 46  0 

455  CONTINUE .. 

460  RR*T*0* 

DDEN*SC<MM,1.MTYPE)-8C(MM-1,i,hTYPE) 

IF (DDgN)464,471, 46  4_ 

464  RR*T«(TLi'-SC(MM-l,l,OTYPE))/DDEN 

471  SCAPt  N 1 *SC ( MH-1 ,2,NTXPE)lRBAT*(SC{8Hl2»MlYPEI“SC(HNilt 

i?;mtypej.>  

76  If  (MTAG(M>-2)  8o. «0*81 
80  R*TI0*EE(2) 

...  ..  GO. TO  82  . 

61  RATI0*6E(2)#EE<3)/EEU> 

82  XX*EE(1)/EE(3) 

YY#1j . 

IF  (MTAG(N)-l)  83,83*04 

83  XX*YY 

84  IF  (HIAGtN)-Jl  86^66*85 ... 

85  YY.XX 

86  CONTINUE 

UU*YY-EEl21*SATl0 . 

VV«XX-EE<2)*RATI0 

UV*EE(2)«(1,*RATI0> 


00003210 

-000«32ti 

00003212 

OOOo322o 

000032IO 

000032*0 

000032(0 

0*<O33|fl 

000o33to 

000o33(o 

OC0q33|o 

000033f 0 
00003350 
OQQo33|o 
0C00337q 
00003360 
OCOo33fo 
OC0034*0 
OCSo34to 
C00034IO 
0000343c 
ocoo34*o 
OC003450 
OCOo34Bo 
OCOo34To 
)OC0034#0 
OC0034»0 
00003500 
OC00351c 
OC0035!o 
0CO0358c 
00003600 

OCOo355o 

OCOo35§o 

OCBo357p 

000035*0 
OCOo35f o 
0000360 0 
OCOo36f o 
00003620 
OCQq363o 
000o36§o 
OCOo365o 
OCQ03680 
OCOo36Tc 
OC0o386o 

000036*0 

0C0037*0 

OCOo37Io 

000037*0 

0000373c 

0C0037*0 


BIO 


3U8T  Ox  08-11-78  15.61Z  

CotJH»EE(i)/<Vv»UU-UVM2> 

C ( 1 • 1 ! Rll)»COMM  _ _ . _ 

C(1.2)*UV»C0HM 

C<1,3)S0. 

C(9,1)RC<1,2) 

C( Si  2 ) *VV»COMM 
C(f.3)«0.' 

- - C(Sd)&C. 

C(S.2)r0. 

C(3i3)sEE(1)/(EE( 1 >/6E  <3)*i.*2.«EE<2J) 

THST*«SlG(N,3)/57.296 

SSaSlNtTHETA) 

CO«COS£THETA) 

S2iSS»SS 

C2»CO*C0 

sco=ss»co 
00  87  1 ] = 1,3 
Do  87  Jjsli3 

87  ref i.jj)=cui,jj) 

55(1,1)802 — - 

D(i.2>=«2 
n<i.3)*sco 
0(2,1)  PS? 

D(2.2)SC2 

D ( Si  3 ) S-SCO  ' 

- D(l=li3r_2,*SCQ > - 

D( 3 • 2 ) R-D (3,1) 

0(3,3) =C?-S2 
DO  88  11=1.3 
DO  88  JJ  = l , 3 
H( ill JJ)=0.0 

DQ  88  KK«1,3  . ....  

88  H ( 1 1 1 J J ) sH  ( I I • JJ ) ♦CUI,Kk)*D(KK,JJ> 

DO  89  11=1,3 

Do  89  JJ=1  i3  _ . 

C(!I|JJ)=0.0 
DO  89  KK  = 1 , 3 

89  C(fliJJJ=C(lI,^J)*Dl!«Ki.l.l)*H(KK<Jj) 

DO  100  I Id, 10 

P ( 1 I ) *0 , 0 
DO. ISO  JJ=1.10 
100  S(TI,JJ>=0.0 
Do  lSO  11=1,3 

oo.itii  wjd-tin 

150  STI 1 1 . JJ)  = 0.0 
VOLsO.O 

C*****F0RM  STIFFNESS  F0R  A CST  TBlAtjGLLA*  ELEMENT*** •••««•#••• 

VOL=O.C 

IF(I*(N,3).NE.IX(N,4>)G0  TO  700 

call  stftriu,2,3)  . ..  . 

xc»(R( 1 ) »R  < J)*P(K) 1/3. 

YC»(Z< I)*Z< J>*Z(K) )/3. 


00808770 

OC 008700 

00fio37f0 

OOOO^BIo 

OOflo3B(o 

JJlOoSBIo 

0C003elo 
OC  003#lo 
OGOqSBSq 
OC  00JB80 
OQQo387o 

osoo^9ta 

OQOo389o 
000039fo 
00003VIO 
0000398o 
008o398o 
. _ _ 0Q0o39|o 

OOOo8V5o 

OCBo396o 

OC005’70 

OCOc398o 

00003970 

0000*080 

0080*010 

000040?0 

OC00405o 

0000*0»0 

OC00405o 

OC00408o 

0000*07o 

OCO0*OBo 

000040fo 

0CO0*l#0 

ocoo4ijto 

oooo^ieo 

OC0041*0 
000041»0 
0 C 0 0 4 3-f  0 
000041#0 
oceo^iJo 
CC0041BO 
ocoo4i9o 

OC004lf5 
00004lf6 
OC004Z»0 
0C0042f 1 
OCOo4Z82 
OC00428* 
0CO04Z85 
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\ 

c*****form  stiffness  for  a cst  quadrilateral  element**^**##.** 

GQ  TO  701  

700  XC«(R< I )*R(J>+R(K>+RIL))'4< 

YC*(7(I)*Z<  J>*Z<K>*Z*LJV4. 

RRR ( 5 ) =XC 

ZZ7(S)»YC 

call  sTftri <4,i,5> 

CALL  STrrRl(l,2 , 5J 

call  stftri<2.3,5) 
call  stftri (3.4.5) 

c*'****calculme  unbalanced  LOADS  DUE  to  TEMPERATURE  CHANGE**#. 
C*****and  stress  relaxation****'**.. «*•*««•*•*••*•*••••••••*«*• 

701  TEMP4(TEmPE(NT-TOL0(N) >*EEC5) 

JF(T1ML  _l7Jl*lAflil2jL 

160  TEMP*0 . 

170  continue 

DSIGIl)  *SLG(N.  l>-*C?*§-l.Q(lii?  >*S2 __  - 

DsiG<2)sSlG(N,l)*s2*SIG(N,2)«C2 

DStG(3l*(SIG<N.l)-Sl5(N,2))*sC0 

... D0.1?Q  JJ  = 1j3 

190  DSIG(JU)=-DSIG( JJ>*(C( JJ*1)*C( jj*2) )*TEMP 
Do  200  11=1.10 
D0.2QO  JJ=l»w 

200  P<  1 1 >*P<  I I )*DSIG<  JJ)*ST<  JJ..U  )*VOL 
C*****ADO  SHEAR  stiffness  OF  FQuNDaTICN*»*»*«*******#«*«**T*** 

IF  UX  (N.3)..EQ»1X  UJ  tA)  ISO  JO  510 

C0MM»v0L*EE<4  ) 

S<9.9)=S(9,9)*C0MH 


Do  500  K«l«2 

_ iHRlOdl 

ID*IH*1 
DO  500  1=1. IH 

saDii(*saDiU/£iin*:£i 

p ( n «P  f J ) -P  (I  D ) *S  ( I . ID  > /S  < I Dt  I D ) 

Do.  500  Jsl.  IH  , , 

500  S(J.n  IS  >_-Sj  J 1 1 p US.  (J  Jxl ) 

c*****calculatf  loads  due  to  gravity******* 

5io  CONTINUE 

. IF  (N-NUMOLl  .580.580. 540  - .... 

340  lF('NNAL,GT.l)GO  TO  5^0 

lF(NNAL.EQ.l.AND.NFLAG.GTr2)00  TO  580 

55.0  LMULi.4 

IF( I X ( N , 3 ) .EQ. IX<N.4> ) IMUL*3 

amul*imul 

..  -D0..56Q  . I-aliJUUL 

560  P(3*I )»P(2*I )-PO(MTyPE)*VOL/AMUL 
380  CONTINUE 

_ 110  RETUM  __ 

End 

subroutine  oned 


00004206 

00004300 

0000438} 

000043*2 

OOOO430S 

00004304 

OCOO4305 

0C!q43|6 

00004307 

0000430a 

00004389 

00004300 

OCOO430J. 

0.C0q43|o 

OCOo438o 

000c43*o 

00004350 

00004300 

00004370 

0CBo436o 

0C0043*c 

00004400 

00004*00 

OC0044?o 

00004403 

0000440; 

0000448a 

ocoo44*a 

ocoo44*a 

00804435 

COOO4*0C 

00004400 

0000446o 

0C0o44*0 

00004500 

00804500 


00004383 

0C0045B5 

000045*0 

00004550 

00004501 

0C0i45g2 

00004533 

OC0O45f3 

OCOo4t>Bo 

00804500 

00004580 

ocan45fa 

OCOfl460o 

00804500 


3U8T  01  08rn-7fl  15. ill 

PARAMETER  HXN*576.Mx8*47S  0C0046J2 

_ CHARACTER  fJX7«6(71_._N*7#6i7>  __  _ 00004*S5 

COMMON  MmNP,NUMEL»NUMPCtN/VOU.TEMP.MTTP|, CAT. NUmNiMANAL.NP* OVER  OC0O4*tO 

l,NnT.NCOUNT.TI,DT,DTr,Tl*Tl.,xC.VC,ST<3.lO).TlKUAC27).NUME<i!7),Tlti  .OCOO^PBO 

2TTT ( 15 1,NUH0L.TImE<NNAL.RL8Y* ISC <EX.EY> DASH, aNALI 200 ).NANALTiTuAS TOC Oo469o 


100 


3,RRR!5),7ZZ<5>,NFLAG»TL1 

COMMON  /MATARG7  E <•  30 , * , 8 J ,.RO < 8 ) ,EE (9 ) , MED < 15 ) iC  IC  (20 , A , 0 > , 

lCCl^liNC^etATiSClMiixlUNSC . _ . 

COMMON  /EUEARG/  ] X (MXE75 ) , MT AG( HxE >;S 1G( MXEf 5 ) , TOLD (MXE > . 
1DEIKMXE)  ,DEl2<MxE>  ,DE2l(MXE>  ,Dfc22CHXE>  ,DSIG|6> ,CC0<4> ,CCC<4) 
2,EGS(MXE,2).EPS(0).Sl-'APCMxe)ANE7(7>.NX7,NY7,TEMPE(NXE> 

COMMON  /PRSARG/  IBc(100I.JBC<100J»P*«10Q> 

COMMON  /ORDARG/  R < MXN> , Z(  MXN) . C,R»  MXN)  »UZ  (MXn)  iCODE(HXN),  T < MgN  > 

COMMON  /baNARG/jiDaNDjNVIISLK. B (.129 1^1 130. 00  > 

COMMON  /US4AHG/  I»J-K.S<10.1C),C(3,3),DI3.3>.H<3,3),P(10>.LM(*>, 

1 F C 3 . 3 1 

CO  100  1*1.8  . 

p<  n*o-'o 
Do  100  *1*1  • 8 

S(i,J)80.Q 

MTYPE= I X ( N > 5 ) 

I»iX(Nil) 

J*IX(Ni2)  ...  ...  

DX«R( J>-R(  I ) 

Dy«ZU)-?(  I > 

.JCL»SOR.T.lDX**2»DY»t?J 

C0SA«DX/XL 
S I NA*DY/xL 

CQtfM*EU,2«MIYPE>t£llx5..NIY£E>/Jll.  

S(i,llRCOSA»COSA«COM0 

S(i,2)rC0SA»SINA*CQMB 

Sii^sj^-siiiU. - . 

S ( 1 , 4 ) s-s ( 1 , 2 ) 

S<2,1)*S(1.2J 


OC0O4699 

00«o46Bo 
0C0fl46|o 
0C0046T0 
OCOo468o 
OC0046fO 
OCOo47|o 
OCOo47fo 
OCOflA7|o 
OCOO4  73o 
OC0047*o 

CC004750 
0C0047|0 
OC0047?o 
0C0C47lo 
0C00479o 
CC  00  48lo 
OC004810 

OCOfl48*o 
0C004BI0 
0C0048Io 
OC  00485o 
OC  0O48(o 

0C0fl48T0 
0000488o 
0C0048»o 
0C0o49«o 
0C8049t0 
0 C 00  49Jq 


. S(2,2)SSINA.»S4NA*CPMM. 

OCBo49a0 

OOio49*o 

, 00004900 

....  . ._  OC  00  49|q 

S (2  ■ 3 ! *-S( 1. 2i 

S ( 2 , 4 ) g»S<  2, 2 J 

S (3 . 1 ) nS ( 1 . 3 ) 

S(3,2 J*S(2,3) 

S (3, 3 ) *S ( 1# 1 ) 

Sf3.4)*S(l,2) 

OC0049To 

OC00498o 

. ._ ...  0C0049fo 

S(4,1)SS(1.4) 

S<4,2J*S(2,4) 

S ( 4 < 3 ) *S (3 . 4 ) 

OC  005Olo 
OOOo’Olo 
OCOflSogB 

S(4,4)sS<2.2> 

EP»S!G<N.'1>/EU,2.Mt*PE) 

DxsDX  »£F 

OCOfl508o 

OCOo50l»o 

OCOo505o 

DyjDY*EP 

P(l)»Sfl7l)*DX*S(l»2>*DY 
P(2)*Sf2.il«DX*SL2ii»4llY  . 

OC00508o 

00«o50To 

OCOflSOBo 

P(3)«-P(l) 

P(4I«-P<7) 

oooo5o9o 

OOOQAlio 

B13 


SHOT  01  OOrli-7#  -15 Ju2 — 


RETURN 

End.  

SUBROUTINE  STFTRUlUJJ.KK) 
PARAMETER  MXNa576,MXE«47* 
character  nx7»6<7).nX7«*<7} 


00005110 

00005180 

000o59|o 

001099*5 

0000596c 


STIFFNESS  MATRIX  (S)  AND  T6ANSF0RHAT10N  MATRIX  (B)  FOR  CST  ELEHBA T000o99f o 
AROUHENTS  (1.2.3)  FOR  TRIANGLE  000fl596o 

- FOR  GUADRIUAIERAl  . 0C6o59fO 

OOOO^Olo 

COMMON  NUMMP » NUMEL * Nl,MPC*NN , VOL  » TEBPJMTYBE . L AT  # NUMN , NANAL . NP 1 0 VER  0000*010 

1. NOT.NCOuNT, T| , DI. PTJ. Tl*TL. XC. lC*5Ti3»lfl) . T|MLA127)»NuME(271iT IF , 0000*0t0 

2TTT( 15 i.NUHOL. TIME. NNaI.NlAV* ISCiEX, EV. DASH. aNAL(20oI.NANALT,T|.A(TOOOo60Io 
3.RRR<5).ZZZ(5>.NfOG*TL1  0000*0*0 

COMMON  /ELEARG/  I X ( PXE . I ) .NT AG( MxE) ISIQ< HX6»9 ) iTOLDtMxEI . 0060*050 

lDEil(MXE),DEl2(MxE)iPE2l(MXE). Dfc22(HXE  ) .DSIGX6  > »CC0(4 ). CCC(4 ) 0000*055 

2, EPS(MXE.2),EPS<6).SCAP(MxEJfNE7(7>.NX7.NY7,TEMPE<MxE>  0000605* 

COMMON  /ORDARG/  R(MXN),2(MXNJ,l.R<MXN).UZ<MXN)|CODElMXN>«T(M)IM  OOO06O60 

COMMON/LSAARG/  INt jN.KNiS (10.10 ) »C( 3.3) 1DUMXI6) jP (lb) »l.H( 4) . 00006000 

1Dllie^sioN  RR(4).77(4>.D(6.0T:'B(3.i0)  I ANGLE!  4 ( F13  lOl  0C9O*0J5 

• 0000*1*0 

1.  INITIALIZATION  OCOo*llo 

. . 0000*120 

,H<l>*lI  0000*1>0 

,M<2)«JJ  0080*1*0 

-M(3>«KK .0510*150 

0000*1*0 

(R(1)«RRR<1I)  OCOo*17o 


1.  INITIALIZATION 

LM('l>*II 

LM(2I»JJ 

LM(3>«KK 

RR ( 1 ) «RRR  < I I ) 

...  ..  RR(2>*RRRUJ) 

RR(3)*RRR(KK> 

z2(i>«nz(ii> 

ZZ(2)»»Z<JJ) 

ZZ ( 3 ) >ZZZ  ( KK ) 

DO  3°  .1*1  >3 

DO  20  J«l,10 

r t f . J>  ■ 0.0 

20  a(  i.j)  L0.ji 

DO  30  J«l.* 

30  D(i,J)  ■ 0.0 


OF  TRIANGLE 


XI  « COMM/2 , 0 
V0L»V0L*XI 


FORM  COEFFICIENT-DISPLACEMENT  MATRIX  IB> 


DU,3)A(2Z(35-ZZ(1>)/C0MM 

D(I,5)»<zZ<1>-ZZ<2>)7C0MM 


O0D0*l|O 

0000*1*0 

000062(0 

0020*210 

OC90*22o 

0000*200 

0C0j)62|o 

0000*2*0 

000o626o 

0090*2?0 

0000*26o 

0000*2*0 

ocse69(o 

0000*3*0 

0000*300 

!>_>_  000j)*3(fl 
0000*3*0 
0000*3*0 
0090*3*0 
0000*370 
0000*3*0 

0 000*3*0 

0000*4*0 

00*0*420 
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D(2 . 2 ) 5 (RR ( 35-RR ( 2 > )/C0MH 

P ( a « * n l RR  ( 1 > 

D(?,8)«(RR(2)-RR(1»)/CQMH 

D(3,1)*D<2.2> 

D<3.2)*P<1.1> 

D(  3. 3 ) |?E  ( 2*  4 ) 

D(J.4)*D(1.3> 

P ( 3i  (itP?  1 

D(3.6)5D(1»5) 

FORM  OB) 

00  50  1«1.3 

DO  31L  Jil. 3 

K a 2*LH(J) 

L a 2*J 

B(l.K-i)*Da*L:l) 

50  8(!,K)s  D(l.L) 


DO  118  M=1.4 

MMalX(NN.M)  

318  ANQCE (M ) aCODE t HM)/57«2958 
Ur2 

irnX(NN.3>.EQ.Ht<BNa4>)  LL«3 

Do  125  Jal.LL 
l»i.M(J> 

tr  l ANGLE  lL) — t(»E.  IU01.  D0_TD-J.25__ 
S IN AaSlNt  ANGLE ( l ) ) 
C08A«C0S<aN0LE(1)> 


DO  124  K*l,3 
TER  ■ 8(k. IJ-1) 

Bi  «*!J:U.a.T£H*CQSA*aiKxIjJJ 'SANA. 

324  B(k.lJ)  a-TEM*SINA*EXK,  U)*C0SA 

325  CONTINUE 

4’  F ORMlLEHENT ’STlF^lifsB Tiatrix  OB)  TMC)  * <¥> 

DO  l30_Jal.lJ5 _ 

DO  130  K*1.3 

!F(B<K«J) >EQ.  0.0)  GO  TO  330 

00,123—1*1.3 

329  F< I,J)SF( I.J)*C< I,K)*B(KiJ>*XI 
130  CONTINUE 


DO  140  1*1.10 
DO  140  K*1.3 

I F ( BIX*  1 ) . EO.  0 . 0 ) _GQ_TA__UU) 

Do  139  J*1.10 

139  S(t.J)5S(I.J)*B(K,l)#F( Ki J ) 


OC0064IQ 
— OOB064|o 

00006400 

00006400 

OO0O6*BO 

00l0«47o 

. saian 

00006580 

00006500 

00006520 

00006500 

00006500 

00006500 

000065*0 

00006570 

0000650c 

00006500 

00006600 

OC0O*0Ifl 

00006600 

OO0q660o 


0000*600 

00006600 

MflO*6It 

00006600 

00006600 

JUOO670O 

000067*0 

00006700 

— .JUtftftZI*. 
00006700 
00006200 
JJflOfl670o 
000067*0 
00006700 

0000*7(0 

00006600 
0000681 o 
00006800 
00006800 
00006600 

20066050 

0000*660 

0000*6*0 

0000*600 

0000*600 

0000*600 

0200*910 

0000*500 

0000*680 


311 BT  fll.  00:01-70 15.6l7 


140  CONTINUE 

91  TRANSFER  OR  ACCUMULATi  tB>  Ik  (ST)  NAl RIX 

DO  200  J*lj3 _ 

DO  200  J*l.iO 

ZOO  STC!.J)*8T(I,J)*t<l,W> 


~CMScK'~rOR  LAST  CALU  fO*  OUaDRILaTER Al 

IFf  1 1 v6«'.3>  AO  TO  290 

RETURN 

290  DO  2*0  1*1.3 

DO  200  J*l.l> 

260  ST(I.J)  * ST ( I . J)/4.0 

RETURN « 


End 

subroutine  stiff 


PARAMETER  MXN»576.MXB*479 
character  nx7*6<7> .NY7»6<7> 


ooooO’lo 

-ooiftotfa 

OOOoOOlo 

O000O9T0 

OCOoOVIo 

000009*0 

00*070|0 

OOloTOio 

00007020 

00007010 

000070*0 

000070*0 

00007000 

Ml07Jlfl_ 

000o70*o 

000070*0 

00007110 

OC0071(0 

000071*0 

0Cfl07l|o 

000o7l*2 

00007189 


COMMON  NyMNPfNUMELiNyMPC+Nj^OU-TlMfiMTTPE. LAY. NUMN.NANAL.NP, OVER  OCflo7l*c 
l.NbT.NCOuNT,Tr,DT.DTT7Tl*TU,XC.YC;6T<3*lO>.T|Ml.AH27).NUME<27),Tl»-IOOBo719o 
2TTT<19),NUMOI.TIME.NNAL.NL»Y*ISC#EX.EY.DASH,ANAU200>.NANALT*TLAST00007100 


3 . RRR  < 9 ) • ZZZ  ( 9 ) . NPL AG*TU 

COMHoTTTM  AT AR57~  E r3B,6,B),R0(8),EE(9).HED( 15) . C IC1 20.4.0). 

1CC( 4«  8) »NCREEP<  8) > SC130 .2.9) .NSC 

COMMON  /ELEARG/  I X (MxE , * ) i.HTAG<MxE )ISlQ(HXE?9 ) . TOtptMXE ) , 

TDElrfMXE) . DElSt  MxE  f,  *EEH  RXE  i7B*22  ( MxTJTDS 1016  >,CCC<4>  .CCC  ( 4 ) 

2.EES(MXET2 ) .EPS (6) »SCAP (NxBlfNE7  <7> .NX7.NY7. TEMPO (BXt > 

COMMON  /PR8ARG/  IBCUOOI*  JBC(H0l*PB(10j>  _ 

COMMON  /ORDARG/  R(m8n).I(MXN)  .UR(MXN)  .UJt(MXN)  i'<rbD67NXN)*T(MXN ) 
COMMON  /BANARG/  mBa*D. NUMBLK , B< 120 > 1 A< 120 . 60 > 

COMMON  /L84ARG/  I . J.K, S(  10. 10 ) . C <3. 3). 0(3 . 3 J . H ( 3 i_3li P t 10).  15(4), 

1 F?3»3l 

C****» INIT I AUlZATlON* *••*•*•***••• AM*******P*t***8#*#*#P***** 

____R1W1NB_2 

N8*30 
ND*2*NB 

ND2jt2«N£ 

STOP*0v0 
numbLkoo 

DO  90  N«i.ND2 


B ( N > • 0 » 0 
DO  50  M»1,ND 
50  JUN.MjJiL.0- 


c***.*form  stifness  matrix  in  Blocks* 
60  numblk*numblk»i 

NH»NB*tNuMBLK*l> 


****** 


NM*NH-NB 
NL*NM-NB*1 


001 

000071*0 
00007180 
OC0071*0 
000072*0 
000072|0 

OO0072|q 

OCflo72Bo 
000072|0 
00007250 
008fl72io 
000072*1 

OOOo7?Io 

00007200 
00007z*0 
OCOo73flo 
000073*o 
0000732o 
_ 0C0o73Qo 
00fio73|o 
000o735o 
OCOo73Bo 

000073(1 

000o73*o 

000o73fo 

000073»o 

000074»0 


B16 


lilfiT  Oi-  48=11^28  lS.ftiZ. 


KSM!ETp2*NL-2 
Mill 


DO  220  M.l.LAY 
Tl*TIM:T IHLA(M) 

TLi»TlHE-TIHLA<H) 

N2»NUME(H> 

DO  210  N*N1.N2 

.If-J-UtNVfn  210.2lfl»iV 


man 


65  DO  80  1>1,4 

IF  ClXCS'.l  )-NL>  80<7fi.70 
70  IF  (IXlK.D-NH)  90*96, 8g 

80  continue 

60  TO  210 

98  ir.,UxtN;a)-tx<Ni2>>..92A.fl4t2_ 

91  CALL  ONEd 
GO  TO  . 165 

92  call.  qu*d  __ 

:r  (VOL)  164.16471ft 

164  WRITE  (6 ,'20031  N 

C»»«»»ADP  EL6HENT  STIFFNESS  TO  TOTAL  8Tlf  rMMH*«.(»(— »*«♦«»« 

165  I X I N . 5 ) *••  I X < N .5) 

DO  166  1*1.4 

166  LH1I)  **•  I XltLi  JJL'2 

DO  200  1*1*4 
DO  200  K*l*2 

J 1»LH  (.U  *K.rK8ri.m 


ooeo7«io 

000o74»o 

OO0oT4fi 

00007400 

00187410 

iiiztfi. 


00007480 
000o74f0 
01 

■"01 


888881! 
OC0075fo 
00flQ75»p 
)075|J 


000o75lo 
00007550 
OC007>80 
OC0075To 
ooaoT’io 

00Bq75$i 


KK*2«lr2*K 
B<t!»*B(Il)*P{KK) 

DO  20.0  J*li4. 

00  200  L*1 *2 
JJ»LM( J)*L*I 1 ♦1"KSMIET 
LL(2» Jr2»L 


0 C 0 0 75f  o 
OCOo76|o 
0C8o761o 
OO0oT6|o 
000076l0 
_Q.COc76Io 


IF(JJ)  200.200.175 
175  IF(ND-JJ)  180,195.195 

ISO  MRiTE  .t0i2884l_J4 ....... 

stop»i:o 

GO  TO  >10 

195  A(fl,JJ)*A<li;jJ)»8<HK.Ll> 

200  CONTINUE 
210  CONTINUE 

NltN2*l _ 

IFINI-NlmEL)  220.22o.225 
220  CONTINUE 

C*»..*ADD  CONCENTRATED  FORCES  OIThin  8LOCK(...»...»c«.(«.t... 
225  Do  255  N.NL.NM 
K.J.N-KSHIFT 

B(K)*9{5)*gZtN.).._T 

B(K-l)SB(K-l)*UR(N> 

IF  (NCOLNT-NP)  255.250.2f0 

250  IF  .(N-NtMN) .252*252,255 

252  UZ < N) *0 , 

UR ( N ) * 0 , 


00007650 
00007680 
.J.0fio767a 
00807680 
O08076f 0 
0C0O7.71O- 

oc8e77lo 

00007780 
000d77Io 
000077f0 
000077*0 
0090771a. 
000o77fo 
00007760 
OCpo77|0 
00007880 
0C007880 
00007811^ 
OCOO  782 0 
00007880 
OC00769o 

0 000  7850 
OCOo7BBo 
000078|0 
000078*0 
00007890 


3U*T  ai.-.A*sli«U- 
255  CONTINUE 

C«*«*.BOUNDaKY  CONDJT | 0N8»B 
i pRCssuRe  67  C« 
ir  ( NuMPc ) 260, 310,260 

260  DO  300  L«1.NUMPC 

liTBCft) 

J.JBC(l) 


giBinai. 


SZ¥T2(T>»ZlJ))«rP 
DR«(R<J)-R< I) >*PP 

26*  lijl'irKSHirr 

JJ«2*Jri<BHlft 
IP  (II)  2*0. 260. 269 

2*9  |P_  1 1 IrN5>  270.270»2»0 

270  SlNAa0 • 0 
CO8Aal>0 

IFO-CODicIU  ^^.272*272  

271  S]NAiSlN(CODE(l>/9)73) 

CO«A«COS<CODE( 11/97.3) 

272  B( 1 I»1 >aB( 1 1~1 )*(C0S6*DZtStBAaDP) 

BdllmBni  >-(SINA*PZ"COSA*DRj 

280  IF  <JJ»  300.300.269 

_2*5  IP  (JJ-ND)  29Q.29Q.3P0  . 

290  S!NA*0.0 
COtAPl. 0 

I r . ( CQDE(J))  291.2*2*292 

:<J>/ 


291  SfNAsSfS  < CODE ( J 1/5773 ) 

C09 A a COS (CODE (J)/57,3) 

292  B( JJ-l) «B ( JJ-l ) ♦( C0S6»DZ7S INA»DR ) . 

B(  Jj) a®  fjj  > - < S 1 Na#PZtCDSJ*DR  > 

3oo  continue 

Co*w»2. displacement  b.  c. 

3lb  DO  400  MsNL.NH 

IT  <M-NUMN>  315.315*400 

_ _?1*  U*«RI») _ 

N,»«H-1-KSM!FT 
IP  (eoOE(H>>  390.400*316 

316  IF  (CqOE(H>-1.)  31?. 070,(17 . 

317  ir  (CODE (M ) -2 . ) 31»,*9C73l$ 

318  IF  (CoDE(M)-3.»  390,380,390 

37 o CALL  MOD IFY< a;B.ND2,MBAND»N»U1 

GO  TO  400 

3so  call  H0D|FV(A.B.ND2,PBANP.B.U) 

390  U«UZtM j 

N«N*1 

CALL  MODlFV( A7B»ND2,0BANP»N,U) 

400  CONTINUE.  _ . - - 

e«aa*»URtTE  BLOCK  OF  EQUATIONS  ON  DISK*  SHIFT  UP  LONER  BLoCk## 
write  02 ) (B<N),(A(N*M) , M*l( mBAKD) ,N*lfND> 

Do  420.Nal,ND ...  ... 

K*N*ND 
B(N) RB  CK  > 


00067910 

_0096.79Ji 

000079(3 

006o79«o 

OO0(792g 

WOo79(o 

00(079*0 

OOB(79Bo 

0CO079|o 

00007’70 
000(79(0 
000(79(0 
0000*000 
000(8010 
0000*02 o 
0060*000 
.0000*010 
0000*000 
00(0*060 

000(6000 

0000*000 
. 00(0*000 
0000*300 
0000*100 
0000*120 
oc6o*iVo 
0000*100 
0000*150 
0000*160 
0000*170 
OO0O810O 
0000*101 
0000*100 
0060*2*0 
0000*210 
0000*200 
0000*260 
000(8200 
0000*2(0 
OC6o*26o 
00(0*270 
0000*200 
0000*200 
. 00(0*3(0 
0000*3X0 
0000*300 
0000*3(0 
0000*38i 
00(0*300 
0000*300 
0000*3(0 
0000*370 


B18 


aim  01  30=11-79  is+Mxi. . 


eoojo.'o 

Ksl.MIL 

A(N,H)iiArK.H) 
<20  A<K.M)ftO.O 


(•••••CHICK  FOR  LOWER  BLOCK. ••^•••••••••••**4»e***< 

ir  (NMrNUMN)  40. 400*400 
<00  CONTINUE 

IF(STQP)  4?fi-SCfl.4*Q 

<90  CALL  EXIT 
900  RETURN 

2903  FORMAT  <?6H0NEGAT|yE  ARE*  ELEHEfcT  NO*  14) 

2904  format  <29moband  wIoTh  exceeds  allowable  |4> 
END 

subroutine  modifyca.b.nebthband.n;u) 

DIMENSION  a<120.60>. B<120> 

DO  290  M.2.MBAND 

K.N-H.l  __  ...... 

ir  f K > 935.235.230 
230  BU04Btfc)-A<K,M)*U 

_ AiKiMlgfl.O 

235  K«N*M-1 

IFiNEQ-K)  250.240, 240 

240  BtKj.BlKj-ALN.Mi.y _ 

A (N.MjiO . 0 

250  CONTINUE  , 


B(N)«U 

RETURN 

ENfl __ 

SUBROUTINE  bansol 

COMMON  /BANARG/  MM,NuMBLKj.B(120>:0(i20» 
HHiAH 

NL. NN.l 
NHsNN.NK 

REWIND  1 

REWIND  2 
NB«0 

GO. TO  150 

C..... REDUCE  EQUATIONS  BY  BLOCKS*.****..*..*** 

c. i , shift  block  of  equations 

loo  NBbNB.1..  . 

DO  125  N.l.NN 

NM. NN.N 

...  BlNltBlML) 

B(NM  )*0 , 0 
DO  125  K*1.MM 

A(N.M)5A(NM,Hj  ... 

125  A<NM.M)«0,0 

READ  NEXT  BLOCK  OF  E9UATICNS  INTO 

ir  cnumblk-nB)  i5o.20Q.ila 
150  READ  (2)  <B<N),<A(N,M),m»1hMMI.B«NL,NH» 
IF  ( NB ) 20u.100.200 


SJ1SUI8 

0000»4fo 

0000*4*0 

00l0«4li 

oote*4to 

0080*410 

OQOa94Bo 

0000*490 

0000*410 

0000*470 

OOBO*4BO 

0000*490 

OC04«5Jo 

0Cflo*51o 

0000*520 

OCOO*5lo 

00fl0*5«0 

0000*590 

J3iJ0*540 

OOOe*5To 

OCOo*59o 

OCOO*59o 

0000*090 

0C«0*6k0 

C400«*to 

0000*0*0 

0000*0*0 

0(00*050 

0090***0 

0000*OTo 

0000*7|0 
. .0000*7*0 

0000*720 
OC00*7l0 

o^ia*7j0 

OCOO*70i 

0000*7(2 

0000*7(0 

0000*750 

0000*770 

0000*700 

0000*790 

0000**00 

0000***0 

0000***0 

oooo**ti 

ocoo*eio 

0000*6*0 

OCOO**5o 


SUIT  01  08=11-7*.  — 15-4^2 

Ca4aaaaaa3.  REDUCE  BLOCK  0^  EQUATIONS 

200  OQ  300  

IF  (A(N,i)>  225,300,22$ 

225  B(N)«BtN)/A(N.l> 

Dq  *75  L*2.NM  

If  (A(N,L))  230.275,230 
230  CbA<N.L)/A(N,1> 



JaO 

no  2$0  Kal.HH 

_ 4»J*l  _ _ _ 

250  A( 1 , j)RA( I , Jj-C«A(N,K) 
8<!>«8Cl)-A(N.U*B<N> 

Aim*  me 

275  CONTINUE 

soo  continue 

Cr*» ••••_« 4.,  WRITE  BLOCK  QF  EQU0T.IQN&  .CM  Ql&K  2S« 
IF  (NUMBLK-NB)  375. 460.375 
375  WRITE  Cl)  (B(N),(A<N*H>,H«2,HN)IN*1,NN> 

aO-.TO.AflO 

C#«»»»BaCK  SUBST  I TUT  I 0N*»»* •••••••••  ••••••#•<!••« 

400  Oo  450  Mal.NN 
N»NN*1~F 
DO  425  K*2«HM 

l»n*k-i 

<25_  B(N.)ifiCN^JLLNxKJjtlik» .... 

nh«n*nN 

B(NH)*B(N> 

_ *50_A(NM,NB).B<N_>_ 

NBaNB-1 

IF  CNB)  475,500,475 

UWM* (i(N) , i aTRTB  ) 7M*‘2,  & h T7K*  I »Nfi V 
backspace  i 


900  K«0 


AeRTMMWS  In  B * OOOo’lfl 

Kafl  0000’1«0 

Dfl.  600  NB«l,NuHBLK ■ OOflflOlfo 

Do  600  Ml.NN  OC0092*0 

NMaN*NN  0C0092*0 

Ka|^«i OCflo9220 

B (K ) »A  CNHi NB ) 000092*0 

return  oooe’zfo 

End osflfl9?fo 

subroutine  creep  oooo*2bo 

PARAMETER  MXNa576,MxEa475  0C0092»5 

COMMON^NUMN^^NUMEL.N5MRC*Ni\if6LTTEM27NTTRE,  LAY",  mjMN7NjTNAL*NP,  OVER  OcSo92Jo 
l,N5T,NCOuNT,TI,DT.DTT,TiiTL,'XC,TC;8T{3.10),TlKLA127).NUME(27l,T|il>,OOOo’2fO 
21 T T < 15 ) , N UMOt. , T I ME.NRAL.N LOT? ISC.EX.EY.PASH, A NALI20o).NANALTaTLAST 0 0 0 o 9 3 0 0 
3,RRR«5!,2ZZ<5),NFLAG*TL1  0C00*3«5 

COMMON  /MaTaRG/  E(3P.6,8),RO(8),EE(9>,HEDU5),CIC<20,4.B),  0C00’3f0 


OCQqOOOi 

O0flo««l(L 

00008070 

0000*8*0 

OOOgBBfo 

0000**50 

ootoMSo 

OMo**|fl 

0000**90 

0000***0 

ocoo**»o 

0000***0 

0000***0 

OOOq8*Oo 

0000*990 

00008909 

0000*0*0 

0000*090 

000o*02o 

OOOo’OBo 

O000*O8l 

000090*0 

OCOO’050 

000090»0 

OCOfl’OTo 

ocao’oeo 

0000*090 

CC0091#0 

0000*1*0 

oooo9xto 

0000*180 

m$f; 

0000*1*0 

OOOfl’lTo 

oooo’lfl 


Da  600  NB«1,NUHBLK 

Do  600  Ml.NN 
NM«N*NN 

_ Kajj«i _ 

600  B(K)*ACNM,NB) 

RETURN 

End 

subroutine  creep 

PARAMETER  MXNa576,MxEa479 


-Character,  Jim* 

COMMON  NUMNP.Nl. 


JUST  01  O0rll-70  15.0lZ 


lCC(4<8)  , NCREEP(O)  eSCt3O,2,0I*NSC  0000*3(0 

C5MM0N  AELEAfiG/  I*l8JUuf.)*8lAG(MXi)l»iaW«E»5J»TQU»iMXiJ* 0010*300 

lDEll<MXE>.DEl2<MxE>.DE2l<MXE>,Dfc22(MNE>i0SlG46»,eCQ<4)«CCC(4J  0000*3*0 

2, ESS (MXE  . 2 > >EPS (6 ) f SCAP (MxE7fNE7 ( 7 > # NX7, NY7, TEMPS (HX0)  0000*300 

COMMON  /PRSARG/  !BC< 100 » . J0C( lOO).PR(lOO)  0000*300 

COMMON  /ORDARG/  R(MXN)iZ(MXN) >LN<MXN) > UZ ( HXN) • CODE (HXN) * T ( HXA ) OO0O*3ro 
COMMON  /BANARG/  hBaMD,NUMBLK(B<120);a<120»60>  0000*300 

eOHNfl^/JJl»RB/  l» 

1F(3,3)  0000*400 

CHARACTER  DASH«6,0VER*6»ETti*l(2)  000o®«*0 

NPRINT^O  ““ 

Nisi 

I Cal 

DOOOO-MalxLAJf 


N2aNUME(M) 

00  550  MM»N1»N2 

n.mm  _ 

TLaTIMrTlMtAIM) 

TLiaTlME.TIMLA(M) 

C»9«>«Emai UAlE  ELEMENT  <TRBSA6S»».««»«».a««aa«aa«a««aahaa.B»m. 

C ki  i cTDcct 

!F  ( IX tN . 2 )- Ix(N , 3 • ) 255*1049255 

255  MTAOUlltl. 

IF  (DStG(4) ) 104,104*299 


259  ir.(0sU(4>«Dsia<5)>  2«ol2to;aoi 

260  MTAQINltg  


GO  TO  104 
261  MTAG(N)<3 

265  »F 

266  MTAGTnIm 

104  iriNFLAG.QT.lJGO  TO  106 
. If  iMJ»Rlm_Afl0ilOljLll0_. 

105  IF(  ISC  i'EO.  0 )GO  TO  750 
HR!TE(6»2005>I1T 

HRiTE(6.2000)LAY.MANAL,TlME.MCOLNT 

WR I TO ( 1, 2007 >LAY, NAN6L1 TIMBTdT 
GO  TO  Z60 

750  MRiTEtS. 2005)07 


000094(0 

0000*400 

0000*400 

Oj)00*4(o_ 

0000*400 

0000*470 

0000*400 

0000*400 

0000*401 

JULA09.400. 

0000*900 

0000*5(0 

0000*9(0 

0000*500 

smisu 


WRITE (6, 2002 >LAY,NANALi TIME, NCOtNT 
760  WRiTE(3)LAY,NANAL.TlBE,NC00NT 

HPXlNTaSJL 

106  CONTINUE 

(•••••calculate  element  strains«9»»«*«»»»*««*9«***«*«*»***0«« 

EP§(115tl./EESm.lI>«ttl.*B0»IMf2T»«5VRBtl0H>-0BSlM^» 

1«(1 .*EES(N»2) >*DSIG(2> I 

EP|(2>4(i,/EES(N.l>)*((l.-BES(N|2>P*a>*DSI6(2>-EES<N*2) 

lEPO<  3)  * ( 2^*U  ♦tlil  N ! i I ) •¥sWj>  )/E«lfN*l) 
EPOTl*tEPS(l)*EPS<2) >/(, 

EP«T2a.9*(OQRT< ( (EPSU)-BP8t2> >«»?>«SP8<3)»»») > 

EPS(4)AEPST1*EPST2 

EP0<5)*EPST1-EPST2 


0000*900 

0000*970 

man 

0000*905 
0.0.00*6(0. 
oeoo*«to 
0000*«(0 
0000*000 
0000*000 
0000*000 
_O0 00*610 


OOOo’OTo 

0000*000 

.0000*090 

0000*700 

0000*701 

0000*7(0 

0000*715 

0000*7(0 

8fUS>l! 

0000*700 

00009700 

0000*700 

oeoo’770 


B21 


lll*T  01  0®Ul'7* 


SDtrf»6PS(l)-EPS(2)  . . . 

EP*  < 6 > « ATAN2  !fP|  < 3 >_,  9fiif  LU1± 

EP8<6)*S7,296»EPS<«> 

DO  700  I IE*1 >5 

7oo  EP4<  i iE>.EP$uiei*joeoeoix 

Ea8REV*EES<N.1>/1O0O0OO. 
iriNFLAG.QT.lTGO  TO  ?5l 

DQ^759  Jflji  „ 

ir«EPSlJUl755.*55i756 
759  ETRCJL>*1H 


GO  TO  750 __ 

756  ETR<JL>*lH* 

759  CONTINUE 

COMPARE  8TRA1M.JO  £T«Al*l_CAP-ACi:iT  IF  APPLICABLE****##** 
!P< 1SC7EO.O)GO  TO  3o9 
lF«EPSC4>-0. >706.706*705 
709  NSX* ( EpS ( 4 ) /SCAP (N  7 ) *100 • *9 . 9 

ifinsxtge, iooojqo  To  720 

IF{ NSX.LT ,79)00  TO  729 
!FCJi»XlfiEJ>aiflQ_TO-Jl6  _ 

ENC0DE»Kx7(IC).729>"  ".NSX."*  • 

729  FoRHATtAl.l3.A2) 

GO. TO  Z 09  . 

726  ENC0DEtNX7<IC>.*29>"  «,NSXi***" 

GO. TO  709 

725  ENC0DE»Nx7<IC1_.72?M  J'lNSAa'L.," 

GO  TO  709 
720  NX7< IC)*OVER 

GO  TO  709  _ . _ . 

706  NX7( JC)*DASH 
NSXaO 

MYtllfljKL— 

Go  TO  *4C 

709  ENC00ECNV7(!C).736>"  ".F.PSU) 

736  FoRMATlAl,F5.1> _ 

740  netUCUn 

IF (1C ,E0 . 7 .OR . N .EQ <NWHE (LAY ) AGO  To  710 

1C*IC*» 

GO  TO  76i 

730  HRtTE(7,742><NE7(I»,Nx7<|),.NY7(I);i*l.:C) 

742  F0RMATt7(lX.l5.2(A*)») _ . J 


ICft 

761  HR!TE(6,2001)N.XC.TC»(DSlGtI)i I *1 J6 )»EPS( 1> « ETR ( 1 ) » EPS! 2 ) • 

. lETR(2)f(EP0( 1 ) > l»3t6> «EABRE>*EEItM»2t«5CAP(N) 

WRITE(S)N.XC.TC, <DSlG<T)il*lf6MEPS<I>.I«l#6). 

1(EGS(N*I). I*1.2).SCAP(N) 

305  WRITE ( 6^2003 )N.XC.TC»( 0 S I GTl) 7 1* 1 1 i )T <EP S (1  >VlVi>6 ) .EABREV 


l.EfS(N*2) 

salts  <: 


LA#JL*1lAL 


C*****RELAX  STRESS  FOR  C«EEP  HjTN.  CONST*NT  STRAIN-FCR  NCBEEPkO 


000 
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0280 
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00 010 3(0 
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l03*o 
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l03*o 
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lOSBo 
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OQOi039o 
OOQiOSBo 

OOOlOSTo 
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3U8T  01  08rll-78  15,41.7 


'•••••••••••••••• »•••••••« 

751  lElUCOyM-HPJiSJLJftilfl 

50  lF(IX<N,2)-lX<N,3>>5V,550,55 
^••••interpolation  of  creep  constants* 
55  I«iX(N«l) 

J* I X ( N»2 > 

K«ix(N*3) 


!X<N,5)»tA9S<IX<N,?>* 

MTVP6«lX(N,5) 

TQLO<N)*TEHPE<N) 

IF  (NCREEP(HTVPE))  290.250>60 
60  NCNaNCREEP(KTYPE) 

! fJl  nil  ?QiA10»12g 

110  TL«TEMPE<N) 

120  DO  1^0  NN*2 > NCR 

IF  (TL-TTT(NN))  125,150.140 
125  TM»TTTtKN)-TTT<NN-i> 

DIFF*  TL-TTT<NN-d 

00  HO  KK»1.4 t : 

130  CCO<KK>*CIC(NN-l,KK,tnrPE)*DirF«<ClClNN,KK,MTTPEJ-ClC<NN 
IE ) >/Tm 

Go  TO  160  _ _ 

140  CONTINUE 
150  DO  155  KK«1.4 

155  CCniKKlaCICtNN.KK.HTXPE) 

160  DO  1*5  KK*1 , 4 

165  CCC(KKI»CC(KK,MTYPE> 

C»****SeLECT  APPROPRIATE  CONST  AN Ts«»»« 

IF  ( DSl G ( 4 ) ) 170.170*175 
170  KK«1 

..  . GO  TOL  U 0 - 

175  KK«3 

180  CC01«CCO(KK) 

CC02SCCO<KK±1)..  - 

CCO3«CC0l 

CC04RCC02 

CCfcllCESUtKl— 

CCC2»CCC(KK*lT 

CCC3»CCCl 

CCC4"CCC2  . 

IF  <DSIG<5>>  l85.lS5«l70 
185  CCO3«CC0(l) 

CC04»CCQ(2I 


CCC3»CCC ( 1 ) 

CCC4»CCC(2> 

. 190  CONTINUE.  „ . ..  . 

^••••MODIFICATION  OF  STRESSES  to  ALLOW  for  CREEP  RELAXATION#* 
C***»*OF  STRESS  AT  constant  STRAIN  on  TNE  APPLICATION  OF  Af ••• 
C*»»»*T ] me  I NCREMENT ••••♦•••••>••  •••••••F*>««»8*»**»»*W»f»* 

THETA* (ESIG(6)-SIG<N*3> 1/57.296 
CO»COStTHETA) 


000i03F2 

ODilOJiO 

OOOlOJfo 

000i03*i 

0C010480 

006l04Ko 

000104CO 
. J0fi0l04|o 

OCOiOARo 

00010450 

00010470 

00010400 

008i04»o 

OOfilOSSo 

00010510 

OOOiOSto 

000l05|o 

OOOlOSRo 

00010500 

_ DQ010560 

-l.KK.NTTPOOOiOSTo 

00010500 
00010500 
000l06»o 
OCO1O6I0 
_000l06*o 

OCO1O680 
OCO1O680 
0CDi06*i 
00010650 
OCOIOASO 
. . 000l067o 

OOO1O680 

000l06»0 

_ . 000l07j0 

OOOl07*o 
000l07*o 

__0Cfll07|0 

0001071b 

008l0750 

000l076o 

000l07T0 

00010780 

000l07*n 

OCO1O8I0 

ocoioeio 

000l08*o 
OOOlOBti 
oooioe*2 
QOOioega 
ocoioeBo 
ocoioeno 
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IllDT  01  06xii-7J 15JHJ 


SS.SlNlTMETA) 

ca.co.co , 

S2.SS.SS 

Siai.ca.<SlG<N,4).Sl*(N,l>>.S2*|Sj«<N.5>.SIG»K»2J> 
S!02.Sa.(SIG<N.4>.SlS(N,l)>.C?.(S]«(N,5).SIG»N|2>) 
ir  (TiT  i99ii*i»i« 

191  DE11 IN J*DS tG<  4 >*CCC2»CC02 

DEi2<N>«DSlQl4)»CCCl«CC01  _ 

DEfl  I NTiDSWyJ  .CCC4«C€6* 

DE*2(NJ.DSIG<5).CCC3*CC0S 

* 192  DE^  11  .il.biil I N > .*2«DE2» ( B ) 

0ET12.  C2»DE12(N >*S2»DE22(b) 

32.DEn(N).C21I!E2lim 

DET22.  S2.DE12(N)*C2*DE22(C0 
DE 1 1 1 N ) » DETli*(DSl0<4>-SlGl).CCC2«CC02 
DI=l2JN).  DETl2.(DSlGU)-S]Gl).CCCl«CCOl 
DE?1 IN).  DET2l.(DSlG^5j-SlB2).CCC4aCeo4 
DE22I N ) « DET22.<DSlGtS)-6lG2)'CCC3>CC03 

SIG(N.l)^ 

SIO(N,2)=0. 

tr«0.  

220  TR.TN.0TT 

DELTilS(DEll(N).OEl2<Nn*OTT 

DElTA2S(hggUN>«PE22lN))»BTT . 

SIQl»D6lTAl*F<l.l) 

SI02«DELTA2.F(2,2) 

Slfi  (_N^1 1 *S  LG.I  N 1 3U  - f Ljl _ 

SIO(N,l)sS!G<Ni2)-Sl52 

C.....CALCULATION  OF  CREEP  RATES  FOR  THE  NEXT  TIME  INTERVAL*.* 

dIi2In1»  Eli2(N)Icit^DTT?cScf > 

DE21<N>«  DE21(N)*(1.-DTT«CCC4> 
DE22fN)«D£22lN)*UA-DTT.CC53J 

oeu  < n J .DEiitN) -s T6i#cccr.cdo2 

DEi2IN>.DEl2(N)-SIgl*CCCl.CC0l 

DE21(N}«DE21<N>-S1.2fCCC4.6C04 

DE22 I N1TDE22 ( N ) -sT52*CCC7iCC03 
C IF«TR-DT)220»250»250 

_250  SlGjN.3>eDSIG<6) _ 

S10<N,4),DSIG<4> 

S!0lN,S)«DSIGI5) 

5*0  CONTINUE 

N1.N2.1 

IF  (Nl-MJMEL)  600. ‘OP. *50 

...600  .CflNllMUL - 

*50  continue 

RETURN 

2005  FqrhaT tlHl PRESENT  time  lNCREHENT..DT-"tr7;2i«.  day3"i  ... 
2000  FORMAT!"  NO.  OF  LIFtS/OR  TIRE  SPAnS-9,15/"  ANALYSIS  ", 
l-CTCLE  NO.".9(1H-).I5,20X.',.KEY  FoR  CHARACTER  POSITjON" 


OOOiOOfo 

OOQiOSBo 

OCOiOOTo 

ocojoeeo 

ooojoafo 

000l09|o 

000109(0 

000l09f 0 
000l09io 
000l09lo 
000l095o 
000i09|o 
000l097o 
000l090o 

000i09fo 

000110*0 

oepiioio 

ooeiioeo 

ooeiioio 

OOOllOSo 

OC011050 

OOOilOEo 

OOOilOTo 

OCCllOOo 

OOflilOfo 

ocoiuio 

OCOillIo 

OCOillEo 

OGOlUBo 

OCOlllRo 

OCOiUfi 

OGSiUDo 

OCOlllBo 

OCOillTo 

ocoiupo 

OCOlUfo 

000il2f0 

0COll2(o 

000ll2to 

OOOllZSo 

0C0ll2*C 

0C0ll2f0 

00011210 

OCOilETo 

0CQ112B0 

OC011290 

OCQxlSfO 

ocaiisio 

OCCllSfo 

000il3|o 

ocoiisio 

OCSilSfo 
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1 


suet  01  00=11-70  is.fiii — 

2,"  FOLLOWING  EPI-X  t EPI«Y  BELOW** /"STRESS  AnC  0TRAINSY. 

3"  Af W-’—'xF IlZxI  MliiliJUl!!  « T1NS1LE  *IRA1H"Z 

4"  APPROXIMATION  NO."«10(1H-J*15,2oX,9BI,ANK  ■ COhpR,  St.IAlN", 
5//"  ElEH". 

6"  CENTROIDtIN. , rilS\lh-^"PSl"4a8il«->*"  AnSL|  "ii.3trMr). 
7"MICR0STRAINS»,14(iHt).«  ANGLE  E<CI  PQIS  S.CP*"/ 

8"  NO.  X-0R3  Y-ORD  SlG-Jl  SlG-Y  SJG-XY  S1G«i  $10". 
...  l!i|  WBT.X  RP1-X  _ 5? ] -X  EP  I-J<Y  *Pi;l  Efll*  JRT.g"* _ 

1"  «10**6  RT10  U- IN  < "// } 

2002  FORMAT!"  NO.  OF  LI^TS/OR  TIME  SPANS-  "*lf/"  ANALYSIS  ", 

1- CYCLE  N0,"1?(1H-),J5/-  ST8ESS  AND  STRAIN  AFTER---", C7’.'2i 

2-  DAYS"/"  APPROXIMATION  N0f«»10 ( 1H- )1 1 5//"  ELEH  ", 
3"C«NTR0ln-IN.  «,l8(lH-)i-PSI"»l0flHi»>."  ANGLE  "»13tlH->)i 
4«MiC»OSTl5AlHS;it4ilH*)A."__.AaGL£  * £iCl  M15S"/"_.Nfi. 

5-X-ORD  Y-ORD  SlG-X  S1G-.Y  SlS-XY  SIG*1  $IG«2  *. 

6"WRT*X  EPI-X  EPI»Y  EPI-XV  EPI»1  EPI-2  WRT*X  *, 

7-*i0"«6  RATIO"//!.. 

2003  F ORMAT il4»lX.2F7.1*3E8,liF7,2|5F0,l#F7.2»F7.3«F5j2> 

2001  F0RMATn4.lX.2F7.i.5P8,liF7.2,F«.l»2<Al»F7.1>.2F8.i, 

lFT’.2.F7.3.rg.g,F7.S) ________  __  

FORMAT!///-'  LIFT  NO.  ",  18/"  LIFT  ANAL  VD76”""* 1 8/ 


ociiisse 
OOAlUfO 
ocoiiaio 
008113*0 
OCOil*0O 
000114*6 
000114(0 
00011406 
000*1406 
00011406 
00011406 
000114*6 
00011406 
000114*6 
00011*00 
OC  011**3 
00011**0 
OCOilSEo 
0001150c 

KftSB 

000il5*o 
00*11*80 
OO0H**e 
000114*6 
00011002 
0001160; 

COMMON  MlMNP,NUMFL*NUMPCiN,.VOL#TEMP.MTyPE,LAY,NUMN,NANAL/<NP.OVER  OCOllOlO 
l,NOT,NC0uNT,TI,DT,DTTJlT^TLl.XC,TC.ST<3.10).TlMLA«27);NvME(27>,Tl(>,OC0ll6f6 
2TTT< 15 ),NUMOL; TIME >NNAL, NLA Y,  jSC.EX.EY, DASH ,aAAL12Oo).NANALT*TUASTOCOi160O 

...  . 00011009 

000110*0 
000110*0 
00011600 
cccii6*0 
00011600 
OCOll0*O 
000117*0 
OCOll7I0 

COMMON  /LS4ARGZ  I / J»K . S < 10 , 1 0 > » C ( 3 , 3 > . D < 3. 3 ) , H ( 3 . 3 >.p UO ) , L» ( * > . OCOll70o 


2007  FORMAT (///- 

1"  TlMElDAYS) ",^10. 2/"  TIME  lNCR(DATS)-».F10,2 

2.7/71"  ELEM  XE(1)  ANGLE.")//) 

END 

subroutine  stress 


PARAMETER  MXN»576.HxE»47> 

CHARACTER  NX7*6<7)»N*7»6<7) 


3.RRR(5>.Z22(5) sNFLAG.TL1 

COMMON  /MATaRG/  E<3?j6,0) iRO{8)iEEl5>rHED(15>,ClC(20,4,e>. 
lCC<4.8).NCREEP(3).Sc!30,2,8T,NSC 

COMMON  /ELEARG/  lX(MxE,S)/.MTAG(MxE)ISIG(MXEr*>iTOLD(MXE>' 
lDEil(MXE),DEl2<MxE,»DE2l(MXE),Dfe22^MXE>»DSIGl6>1CC0(4).CCC<4 ) 

2i EESTmXE , 2 ) .EPS! 6) *SEAP I PXC  ) *N£7 I?T  *NX7 ,NY7 , TEHPE  t BxE ) 

COMMON  /PRSARG/  iSC^lOQ) , J0C < 100 ) -PR ( 100 ) 

COMMON  /ORDARG/  R<mXN1«?(MXN) iLR(MXN) aU? (MXN>iC006iHXN)*T1MSN) 
COMMON  /RANARG/  MB A«D . NUMBLK , B ( 120 > 1 A < 120 . 60 ) 


1 Ff.3i3)  ._.  ......  , „ 

c*«»«*compute  element  stresses* 

DO  50  1*1, A 

50  DSIG(ns,0,0  

IF  (!x!N.3)-IXlN.2>)  90,80/90 
C****»F OR  OnE-D  ELEMENT •*••••*♦*•*«< 
80  I*!X(NU) 

J»IX(N*2) 

DX»R<J)-R(D 

Dy»Z« 

XL«S0R,T(DX*»2*DY»*2) 

Du»BI 2* J-l )-B ( 2*1 -1 ) 


OC0ll7?p 

0C0il78i 

00011700 
OC  011  700 
OC01176O 
0C0il78i 
0C0ll7*c 
OCOll70o 
OC0ll7fo 
OCS1180O 
OCOllBJo 
OCOllSBc 
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SllBT  Oi  _ 08^11-78 13.01? 

DV*B(29j)-B(2*I>  , . 

_ oi..pv«Dm-k_*e^OJ{OJ‘! _ 

DStG<lI«DL*E<1.2.NTyPE>/*L*SlGt8#l>*SIG<!)l,4) 

IF  (NCOUNT-NP)  85.84*84 

84  SI0(N,4)«DS16(1) 

S!8(N,i)*0. 

89  XC»0.0 

YCaO.O  _ 

GO  TB  320 

C*»***F 04  TWO-0  ELEhE NT ••••••••••••••• ••••••*•*•••*••« 

90  CALL  OU*t) _ 

DO  l20  1*1.4 

1 1*2*1 

JJmltllllUil _ 

P( ! !-l)»B( JJ-1> 

120  P(M)*B(JJ) 

DO  150  1*9.10 _ . 

P ( t ) » 0 • 0 
KK.J-1 

DO  1?Q  11*1. KK. 

150  P(t >»PtI)-S(I.K>*P<K* 

D(l.l)*Ci. 

0(2.1)40.  - - 

D(3,l)*0.0 
DO  170  1*1.3 

DQ_170  ... 

170  D<I,1)4D<I.1>*ST<!.K>*P(K) 

ThBTA*  SIG(N.3>/57.2?6 

. . CqaCOSSTHEJAA. 

SS*sIn(ThETA5 

C2*C0*CC 

S2*SS*SS  _ _ . . . 

sco*ss*co 

DSIG<H*  SIG<N,4>*C2*SIG(N,.5)*S2-DSIG(1) 

DS!G<2>*  SIG(.N,4)*S2*SIG(N.5)*C2-0SIG(2) 
DSrG/3)*<S}G(N74>-Sl*rN,T))*SCO*DslG<3> 

Do  ISO  1=1.3 

00.180. K«l»l ..... 

180  DS!GU>*  DSIG<I>*  C< r,K)*D(K.l> 

C****«OUTPUT  STRESSES*. 

c**... calculate  principal  stresses*.******#*********** 

AA*(DSIG(1)*DSIG(2>)^2. 

BBs  (DSJG(1)-DSIG(2) >/2. 

CR»  SQRT(BB**2*DSIQ(3)**2) 

DS I G( 4 ) * aA*CR 
DS!  G (5 1 * AA-CR 

IF  <«8B,eO.0-O>.AND, IDSJgtSj.feo.O.P))  GO  TO  320 
DSl'G<6l*ATAN2<nSlG<3»  .BBl/2.’ 
DSiG(6>*57.2«6*DSlG<*) 

320  RETURN 
END 

8 EXECUTE 
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oooiiero 
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00011920 
000ll9ti 
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OOBll9lo 
OOOll9fo 
OCOil9fo 
OCOil96o 
OC011970 
OCOilVBo 
ocoii’Bo 
0001200 0 
000l20*0 

OCOi20to 

CCOi203o 

0003.20*0 

OC0l2O5o 

000i208o 

OCOl20To 

0C012000 

0C0l20f0 

OCOi210o 

0C012110 

0C012120 

000l21fo 

0C0l2ll0 

OCOi215o 

O00l21Bo 

00012170 

0C0i2lJc 

OC0l21fO 

OCOl22|o 

OCOl22«l 

0C0l2202 

OC012210 

000l22fo 

0C012280 

0CDl22«0 

OC012250 

0C0l226e 

00012270 

OCOl22Po 

0C0l2Z»0 

00012300 

000l23Ic 
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